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Abstract 

The status of experimental tests of general relativity and of theo- 
retical frameworks for analysing them are reviewed. Einstein's equiv- 
alence principle (EEP) is well supported by experiments such as the 
Eotvos experiment, tests of special relativity, and the gravitational 
redshift experiment. Future tests of EEP and of the inverse square 
law will search for new interactions arising from unification or quan- 
tum gravity. Tests of general relativity at the post-Newtonian level 
have reached high precision, including the light deflection, the Shapiro 
time delay, the perihelion advance of Mercury, and the Nordtvedt ef- 
fect in lunar motion. Gravitational wave damping has been detected 
to half a percent using the binary pulsar, and new binary pulsar sys- 
tems may yield further improvements. When direct observation of 
gravitational radiation from astrophysical sources begins, new tests of 
general relativity will be possible. 



1 Introduction 

At the time of the birth of general relativity (GR), experimental confirmation 
was almost a side issue. Einstein did calculate observable effects of general 
relativity, such as the perihelion advance of Mercury, which he knew to be an 
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unsolved problem, and the deflection of light, which was subsequently veri- 
fied, but compared to the inner consistency and elegance of the theory, he re- 
garded such empirical questions as almost peripheral. But today, experimen- 
tal gravitation is a major component of the field, characterized by continuing 
efforts to test the theory's predictions, to search for gravitational imprints 
of high-energy particle interactions, and to detect gravitational waves from 
astronomical sources. 

The modern history of experimental relativity can be divided roughly into 
four periods. Genesis, Hibernation, a Golden Era, and the Quest for Strong 
Gravity. The Genesis (1887-1919) comprises the period of the two great ex- 
periments which were the foundation of relativistic physics — the Michelson- 
Morley experiment and the Eotvos experiment — and the two immediate con- 
firmations of GR — the deflection of light and the perihelion advance of Mer- 
cury. Following this was a period of Hibernation (1920-1960) during which 
theoretical work temporarily outstripped technology and experimental pos- 
sibilities, and, as a consequence the field stagnated and was relegated to the 
backwaters of physics and astronomy. 

But beginning around 1960, astronomical discoveries (quasars, pulsars, 
cosmic background radiation) and new experiments pushed GR to the fore- 
front. Experimental gravitation experienced a Golden Era (1960-1980) dur- 
ing which a systematic, world-wide effort took place to understand the ob- 
servable predictions of GR, to compare and contrast them with the predic- 
tions of alternative theories of gravity, and to perform new experiments to 
test them. The period began with an experiment to confirm the gravita- 
tional frequency shift of light (1960) and ended with the reported decrease 
in the orbital period of the binary pulsar at a rate consistent with the gen- 
eral relativity prediction of gravity-wave energy loss (1979). The results all 
supported GR, and most alternative theories of gravity fell by the wayside 
(for a popular review, see ||144|| ). 

Since 1980, the field has entered what might be termed a Quest for Strong 
Gravity. Many of the remaining interesting weak-field predictions of the 
theory are extremely small and difficult to check, in some cases requiring 
further technological development to bring them into detectable range. The 
sense of a systematic assault on the weak-field predictions of GR has been 
supplanted to some extent by an opportunistic approach in which novel and 
unexpected (and sometimes inexpensive) tests of gravity have arisen from 
new theoretical ideas or experimental techniques, often from unlikely sources. 
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Examples include the use of laser-cooled atom and ion traps to perform ultra- 
precise tests of special relativity; the proposal of a "fifth" force, which led 
to a host of new tests of the weak equivalence principle; and recent ideas of 
large extra dimensions, which have motived new tests of the inverse square 
law of gravity at sub-millimeter scales. Several major ongoing efforts also 
continue, principally the Stanford Gyroscope experiment, known as Gravity 
Probe-B. 

Instead, much of the focus has shifted to experiments which can probe 
the effects of strong gravitational fields. The principal figure of merit that 
distinguishes strong from weak gravity is the quantity e ~ GM/Rc^, where 
G is the Newtonian gravitational constant, M is the characteristic mass scale 
of the phenomenon, R is the characteristic distance scale, and c is the speed 
of light. Near the event horizon of a non-rotating black hole, or for the 
expanding observable universe, e ~ 0.5; for neutron stars, e ~ 0.2. These 
are the regimes of strong gravity. For the solar system e < 10^^; this is 
the regime of weak gravity. At one extreme are the strong gravitational 
fields associated with Planck-scale physics. Will unification of the forces, 
or quantization of gravity at this scale leave observable effects accessible 
by experiment? Dramatically improved tests of the equivalence principle 
or of the inverse square law are being designed, to search for or bound the 
imprinted effects of Planck-scale phenomena. At the other extreme are the 
strong fields associated with compact objects such as black holes or neutron 
stars. Astrophysical observations and gravitational-wave detectors are being 
planned to explore and test GR in the strong-field, highly-dynamical regime 
associated with the formation and dynamics of these objects. 

In this living review, we shall survey the theoretical frameworks for study- 
ing experimental gravitation, summarize the current status of experiments, 
and attempt to chart the future of the subject. We shall not provide complete 
references to early work done in this field but instead will refer the reader to 
the appropriate review articles and monographs, specifically to Theory and 
Experiment in Gravitational Physics ||145|| , hereafter referred to as TEGP. 
Additional recent reviews in this subject are ||139|, |143|, |147|, |149|, BBI, |117 



References to TEGP will be by chapter or section, e.g. "TEGP 8.9" 
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2 Tests of the Foundations of Gravitation 
Theory 



2.1 The Einstein Equivalence Principle 

The principle of equivalence has historically played an important role in the 
development of gravitation theory. Newton regarded this principle as such a 
cornerstone of mechanics that he devoted the opening paragraph of the Prin- 
cipia to it. In 1907, Einstein used the principle as a basic element of general 
relativity. We now regard the principle of equivalence as the foundation, not 
of Newtonian gravity or of GR, but of the broader idea that spacetime is 
curved. 

One elementary equivalence principle is the kind Newton had in mind 
when he stated that the property of a body called "mass" is proportional to 
the "weight", and is known as the weak equivalence principle (WEP). An 
alternative statement of WEP is that the trajectory of a freely falling body 
(one not acted upon by such forces as electromagnetism and too small to be 
affected by tidal gravitational forces) is independent of its internal structure 
and composition. In the simplest case of dropping two different bodies in a 
gravitational field, WEP states that the bodies fall with the same acceleration 
(this is often termed the Universality of Free Fall, or UFF). 

A more powerful and far-reaching equivalence principle is known as the 
Einstein equivalence principle (EEP). It states that 

1. WEP is vahd. 

2. The outcome of any local non-gravitational experiment is independent 
of the velocity of the freely-falling reference frame in which it is per- 
formed. 

3. The outcome of any local non-gravitational experiment is independent 
of where and when in the universe it is performed. 

The second piece of EEP is called local Lorentz invariance (LLI), and the 
third piece is called local position invariance (LPI). 

For example, a measurement of the electric force between two charged 
bodies is a local non-gravitational experiment; a measurement of the gravi- 
tational force between two bodies (Cavendish experiment) is not. 
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The Einstein equivalence principle is the heart and soul of gravitational 
theory, for it is possible to argue convincingly that if EEP is valid, then 
gravitation must be a "curved spacetime" phenomenon, in other words, the 
effects of gravity must be equivalent to the effects of living in a curved space- 
time. As a consequence of this argument, the only theories of gravity that 
can embody EEP are those that satisfy the postulates of "metric theories of 
gravity", which are 

1. Spacetime is endowed with a symmetric metric. 

2. The trajectories of freely falling bodies are geodesies of that metric. 

3. In local freely falling reference frames, the non-gravitational laws of 
physics are those written in the language of special relativity. 

The argument that leads to this conclusion simply notes that, if EEP is 
valid, then in local freely falling frames, the laws governing experiments must 
be independent of the velocity of the frame (local Lorentz invariance), with 
constant values for the various atomic constants (in order to be independent 
of location). The only laws we know of that fulfill this are those that are 
compatible with special relativity, such as Maxwell's equations of electro- 
magnetism. Furthermore, in local freely falling frames, test bodies appear to 
be unaccelerated, in other words they move on straight lines; but such "lo- 
cally straight" lines simply correspond to "geodesies" in a curved spacetime 
(TEGP 2.3). 

General relativity is a metric theory of gravity, but then so are many 
others, including the Brans-Dicke theory. The nonsymmetric gravitation 
theory (NGT) of Moffat is not a metric theory. Neither, in this narrow 
sense, is superstring theory (see Sec. |0| ), which, while based fundamentally 
on a spacetime metric, introduces additional fields (dilatons, moduli) that 
can couple to material stress-energy in a way that can lead to violations, say, 
of WEP. So the notion of curved spacetime is a very general and fundamental 
one, and therefore it is important to test the various aspects of the Einstein 
Equivalence Principle thoroughly. 

A direct test of WEP is the comparison of the acceleration of two laborat- 
ory-sized bodies of different composition in an external gravitational field. If 
the principle were violated, then the accelerations of different bodies would 
differ. The simplest way to quantify such possible violations of WEP in a 
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form suitable for comparison with experiment is to suppose that for a body 
with inertial mass mj, the passive gravitational mass mp is no longer equal to 
mi, so that in a gravitational field g, the acceleration is given by mja = mpg. 
Now the inertial mass of a typical laboratory body is made up of several types 
of mass-energy: rest energy, electromagnetic energy, weak-interaction energy, 
and so on. If one of these forms of energy contributes to mp differently than 
it does to mj, a violation of WEP would result. One could then write 

mp = mj + J2v'^E'^/c\ (1) 

A 

where is the internal energy of the body generated by interaction A, and 
7]^ is a dimensionless parameter that measures the strength of the violation of 
WEP induced by that interaction, and c is the speed of light. A measurement 
or limit on the fractional difference in acceleration between two bodies then 
yields a quantity called the "Eotvos ratio" given by 

_ , / Ji^ _ (2) 

|ai + a2| ^ \mic^ m2C'^ J 

where we drop the subscript I from the inertial masses. Thus, experimental 
limits on i] place limits on the WEP-violation parameters r]^. 

Many high-precision Eotvos-type experiments have been performed, from 
the pendulum experiments of Newton, Bessel and Potter, to the classic 
torsion-balance measurements of Eotvos |S7|, Dicke [ 113 |, Braginsky |^ 



and their collaborators. In the modern torsion-balance experiments, two 
objects of different composition are connected by a rod or placed on a tray 
and suspended in a horizontal orientation by a fine wire. If the gravita- 
tional acceleration of the bodies differs, there will be a torque induced on 
the suspension wire, related to the angle between the wire and the direction 
of the gravitational acceleration g. If the entire apparatus is rotated about 
some direction with angular velocity u, the torque will be modulated with 
period 27i/uj. In the experiments of Eotvos and his collaborators, the wire 
and g were not quite parallel because of the centripetal acceleration on the 
apparatus due to the Earth's rotation; the apparatus was rotated about the 
direction of the wire. In the Dicke and Braginsky experiments, g was that of 
the Sun, and the rotation of the Earth provided the modulation of the torque 
at a period of 24 hr (TEGP 2.4(a)). Beginning in the late 1980s, numerous 
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experiments were carried out primarily to search for a "fifth force" (see Sec. 
|27^ ), but their null results also constituted tests of WEP. In the "free-fall 
Galileo experiment" performed at the University of Colorado, the relative 
free-fall acceleration of two bodies made of uranium and copper was mea- 
sured using a laser interferometric technique. The "Eot-Wash" experiments 
carried out at the University of Washington used a sophisticated torsion 
balance tray to compare the accelerations of various materials toward local 
topography on Earth, movable laboratory masses, the Sun and the galaxy 



121| , ^, and have recently reached levels of 4 x 10^^^. The resulting upper 
limits on t] are summarized in Figure ^ (TEGP 14.1; for a bibliography of 
experiments, see []6l|). 

The second ingredient of EEP, local Lorentz invariance, has been tested 
to high-precision in the "mass anisotropy" experiments: the classic versions 
are the Hughes-Drever experiments, performed in the period 1959-60 inde- 
pendently by Hughes and collaborators at Yale University, and by Drever 
at Glasgow University (TEGP 2.4(b)). Dramatically improved versions were 
carried out during the late 1980s using laser-cooled trapped atom techniques 
(TEGP 14.1). A simple and useful way of interpreting these experiments is 
to suppose that the electromagnetic interactions suffer a slight violation of 
Lorentz invariance, through a change in the speed of electromagnetic radia- 
tion c relative to the limiting speed of material test particles (cq, chosen to 
be unity via a choice of units), in other words, c 7^ 1 (see Sec. |2.2.3|) . Such a 
violation necessarily selects a preferred universal rest frame, presumably that 
of the cosmic background radiation, through which we are moving at about 
300 km/s. Such a Lorentz-non-invariant electromagnetic interaction would 
cause shifts in the energy levels of atoms and nuclei that depend on the ori- 
entation of the quantization axis of the state relative to our universal velocity 
vector, and on the quantum numbers of the state. The presence or absence of 
such energy shifts can be examined by measuring the energy of one such state 
relative to another state that is either unaffected or is affected differently by 
the supposed violation. One way is to look for a shifting of the energy levels 
of states that are ordinarily equally spaced, such as the four J=3/2 ground 
states of the ^Li nucleus in a magnetic field (Drever experiment); another is 
to compare the levels of a complex nucleus with the atomic hyperfine levels 
of a hydrogen maser clock. These experiments have all yielded extremely 
accurate results, quoted as limits on the parameter 6 = — 1 in Figure 
0. Also included for comparison is the corresponding limit obtained from 
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Figure 1: Selected tests of the Weak Equivalence Principle, showing bounds 
on ?7, which measures fractional difference in acceleration of different mate- 
rials or bodies. Free-fall and Eot-Wash experiments originally performed to 
search for fifth force. Blue band shows current bounds on ri for gravitating 
bodies from lunar laser ranging (LURE). 
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Michelson-Morley type experiments (for a review, see |[75|). 

Recent advances in atomic spectroscopy and atomic timekeeping have 
made it possible to test LLI by checking the isotropy of the speed of hght 
using one-way propagation (as opposed to round-trip propagation, as in the 
Michelson-Morley experiment). In one experiment, for example, the relative 
phases of two hydrogen maser clocks at two stations of NASA 's Deep Space 
Tracking Network were compared over five rotations of the Earth by propa- 
gating a light signal one-way along an ultrastable fiberoptic link connecting 
them (see Sec. |2.2.3| ). Although the bounds from these experiments are not 



as tight as those from mass-anisotropy experiments, they probe directly the 
fundamental postulates of special relativity, and thereby of LLI. (TEGP 14.1, 

The principle of local position invariance, the third part of EEP, can be 
tested by the gravitational redshift experiment, the first experimental test 
of gravitation proposed by Einstein. Despite the fact that Einstein regarded 
this as a crucial test of GR, we now realize that it does not distinguish be- 
tween OR and any other metric theory of gravity, instead is a test only of 
EEP. A typical gravitational redshift experiment measures the frequency or 
wavelength shift Z = Au/u = — AA/A between two identical frequency stan- 
dards (clocks) placed at rest at different heights in a static gravitational field. 
If the frequency of a given type of atomic clock is the same when measured 
in a local, momentarily comoving freely falling frame (Lorentz frame), in- 
dependent of the location or velocity of that frame, then the comparison of 
frequencies of two clocks at rest at different locations boils down to a com- 
parison of the velocities of two local Lorentz frames, one at rest with respect 
to one clock at the moment of emission of its signal, the other at rest with 
respect to the other clock at the moment of reception of the signal. The 
frequency shift is then a consequence of the first-order Doppler shift between 
the frames. The structure of the clock plays no role whatsoever. The result 
is a shift 

Z = AU/c^ , (3) 

where AU is the difference in the Newtonian gravitational potential between 
the receiver and the emitter. If LPI is not valid, then it turns out that the 
shift can be written 

Z = {1 + a)AU/c^ , (4) 
where the parameter a may depend upon the nature of the clock whose shift 
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Figure 2: Selected tests of local Lorentz invariance showing bounds on pa- 
rameter S, which measures degree of violation of Lorentz invariance in elec- 
tromagnetism. Michelson-Morley, Joos, and Brillet-Hall experiments test 
isotropy of round-trip speed of light, the latter experiment using laser tech- 
nology. Centrifuge, two-photon absorption (TPA) and JPL experiments test 
isotropy of light speed using one-way x)ropagation. Remaining four experi- 
ments test isotropy of nuclear energy levels. Limits assume speed of Earth 
of 300 km/s relative to the mean rest frame of the universe. 



is being measured (see TEGP 2.4(c) for details). 

The first successful, high-precision redshift measurement was the series of 
Pound- Rebka- Snider experiments of 1960-1965, that measured the frequency 
shift of gamma-ray photons from ^^Fe as they ascended or descended the 
Jefferson Physical Laboratory tower at Harvard University. The high accu- 
racy achieved - one percent - was obtained by making use of the Mossbauer 
effect to produce a narrow resonance line whose shift could be accurately 
determined. Other experiments since 1960 measured the shift of spectral 
lines in the Sun's gravitational field and the change in rate of atomic clocks 
transported aloft on aircraft, rockets and satellites. Figure ^ summarizes the 
important redshift experiments that have been performed since 1960 (TEGP 
2.4(c)). 

The most precise standard redshift test to date was the Vessot-Levine 
rocket experiment that took place in June 1976 ||130||. A hydrogen-maser 



clock was flown on a rocket to an altitude of about 10,000 km and its fre- 
quency compared to a similar clock on the ground. The experiment took 
advantage of the masers' frequency stability by monitoring the frequency 
shift as a function of altitude. A sophisticated data acquisition scheme accu- 
rately eliminated all effects of the first-order Doppler shift due to the rocket's 
motion, while tracking data were used to determine the payload's location 
and the velocity (to evaluate the potential difference At/, and the special rel- 
ativistic time dilation). Analysis of the data yielded a limit |a| < 2 x 10~^. 

A "null" redshift experiment performed in 1978 tested whether the rel- 
ative rates of two different clocks depended upon position. Two hydrogen 
maser clocks and an ensemble of three superconducting-cavity stabilized os- 
cillator (scso) clocks were compared over a 10-day period. During the pe- 
riod of the experiment, the solar potential U/c^ changed sinusoidally with 
a 24-hour period by 3 x 10"^'^ because of the Earth's rotation, and changed 
linearly at 3 x 10^^^ per day because the Earth is 90 degrees from perihelion 
in April. However, analysis of the data revealed no variations of either type 
within experimental errors, leading to a limit on the LPI violation parameter 
— < 2 X 10~^ ||128|]. This bound has been improved using more 



stable frequency standards [^, |109|| . The varying gravitational redshift of 



Earth-bound clocks relative to the highly stable Millisecond Pulsar, caused 
by the Earth's motion in the solar gravitational field around the Earth- Moon 
center of mass (amplitude 4000 km), has been measured to about 10 percent, 
and the redshift of stable oscillator clocks on the Voyager spacecraft caused 
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Figure 3: Selected tests of local position invariance via gravitational redshift 
experiments, showing bounds on a, which measures degree of deviation of 
redshift from the formula Au/u = AU/c^. 
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by Saturn's gravitational field yielded a one percent test. The solar gravita- 
tional redsliift has been tested to about two percent using infrared oxygen 
triplet lines at the limb of the Sun, and to one percent using oscillator clocks 
on the Galileo spacecraft (TEGP 2.4(c) and 14.1(a)). 

Modern advances in navigation using Earth-orbiting atomic clocks and 
accurate time-transfer must routinely take gravitational redshift and time- 
dilation effects into account. For example, the Global Positioning System 
(GPS) provides absolute accuracies of around 15 m (even better in its military 
mode) anywhere on Earth, which corresponds to 50 nanoseconds in time 
accuracy at all times. Yet the difference in rate between satellite and ground 
clocks as a result of special and general relativistic effects is a whopping 39 
microseconds per day (46/is from the gravitational redshift, and — 7/xs from 
time dilation). If these effects were not accurately accounted for, GPS would 
fail to function at its stated accuracy. This represents a welcome practical 
application of GR! (For the role of GR in GPS, see 0; for a popular essay, 
see [Hi.) 

Local position invariance also refers to position in time. If LPI is satisfied, 
the fundamental constants of non-gravitational physics should be constants 
in time. Table shows current bounds on cosmological variations in selected 
dimensionless constants. For discussion and references to early work, see 
TEGP 2.4(c). 

2.2 Theoretical Frameworks for Analyzing EEP 
2.2.1 Schiff's Conjecture 

Because the three parts of the Einstein equivalence principle discussed above 
are so very different in their empirical consequences, it is tempting to regard 
them as independent theoretical principles. On the other hand, any complete 
and self-consistent gravitation theory must possess sufficient mathematical 
machinery to make predictions for the outcomes of experiments that test each 
principle, and because there are limits to the number of ways that gravitation 
can be meshed with the special relativistic laws of physics, one might not 
be surprised if there were theoretical connections between the three sub- 
principles. For instance, the same mathematical formalism that produces 
equations describing the free fall of a hydrogen atom must also produce 
equations that determine the energy levels of hydrogen in a gravitational field, 
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Constant k 



Limit onk/k 
per Hubble time 
yr 



1.2 X 10^° 



Method 



Fine structure constant 
a = /he 



4 X 10-^ 
9 X 10~^ 
6 X 10"^ 
6 X 10~^ 



H-maser vs Hg ion clock [|109|| 
^"^Rb fountain vs Cs clock |TT5| 
Oklo Natural Reactor [|0l 
21-cm vs molecular absorption 

0) 



at Z = 0.7 



Weak interaction constant 
p = Gfmlc/h^ 


1 

0.1 
0.06 


""^^^Re, decay rates 
Oklo Natural Reactor |40 
Big bang nucleosynthesis 


91 


, 112| 


e-p mass ratio 


1 


Mass shift in quasar 
spectra at Z ~ 2 






Proton g-factor (gp) 


10"^ 
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at Z = 0.7 ( 56 


ar absorption 
) 



Table 1: Bounds on cosmo logical variation of fundamental constants of non- 
gravitational physics. For references to earlier work, see TEGP 2.4(c). 
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and thereby the ticking rate of a hydrogen maser clock. Hence a violation 
of EEP in the fundamental machinery of a theory that manifests itself as 
a violation of WEP might also be expected to show up as a violation of 
local position invariance. Around 1960, Schiff conjectured that this kind of 
connection was a necessary feature of any self-consistent theory of gravity. 
More precisely, Schiff's conjecture states that any complete, self- consistent 
theory of gravity that embodies WEP necessarily embodies EEP. In other 
words, the validity of WEP alone guarantees the validity of local Lorentz 
and position invariance, and thereby of EEP. 

If Schiff's conjecture is correct, then Eotvos experiments may be seen 
as the direct empirical foundation for EEP, hence for the interpretation of 
gravity as a curved-spacetime phenomenon. Of course, a rigorous proof of 
such a conjecture is impossible (indeed, some special counter-examples are 
known) yet a number of powerful "plausibility" arguments can be formulated. 

The most general and elegant of these arguments is based upon the as- 
sumption of energy conservation. This assumption allows one to perform 
very simple cyclic gedanken experiments in which the energy at the end of 
the cycle must equal that at the beginning of the cycle. This approach was 



pioneered by Dicke, Nordtvedt and Haugan (see, e.g. fT^). A system in a 
quantum state A decays to state B, emitting a quantum of frequency u. The 
quantum falls a height H in an external gravitational field and is shifted 
to frequency z/', while the system in state B falls with acceleration qb- At 
the bottom, state A is rebuilt out of state B, the quantum of frequency u', 
and the kinetic energy msgBH that state B has gained during its fall. The 
energy left over must be exactly enough, niAgAH, to raise state A to its 
original location. (Here an assumption of local Lorentz invariance permits 
the inertial masses tua and to be identified with the total energies of the 
bodies.) If gA and gs depend on that portion of the internal energy of the 
states that was involved in the quantum transition from AtoB according to 

gA = gil+aEA/rriAC^) , gs = g{l+aEB/mBC^) , Ea-Eb = hi^ (5) 

(violation of WEP), then by conservation of energy, there must be a corre- 
sponding violation of LPI in the frequency shift of the form (to lowest order 
in hu/mc^) 

Z={u'- v)lv' = (1 + a)gElc^ = (1 + a)Af//c2 . (6) 
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Haugan generalized this approach to include violations of LLI ( , TEGP 
2.5). 



2.2.2 The THejj, Formalism 

The first successful attempt to prove Schiff 's conjecture more formally was 
made by Lightman and Lee [Q. They developed a framework called the 
THefi formalism that encompasses all metric theories of gravity and many 
non-metric theories (Box 1). It restricts attention to the behavior of charged 
particles (electromagnetic interactions only) in an external static spherically 
symmetric (SSS) gravitational field, described by a potential U. It charac- 
terizes the motion of the charged particles in the external potential by two 
arbitrary functions T{U) and H{U), and characterizes the response of elec- 
tromagnetic fields to the external potential (gravitationally modified Maxwell 
equations) by two functions e(f/) and /i(t/). The forms of T, H, e and /i vary 
from theory to theory, but every metric theory satisfies 

e = ^^=iH/Ty/\ (7) 

for all U. This consequence follows from the action of electrodynamics with 
a "minimal" or metric coupling: 

- 1^ [ V^9'"'9''^F,^Faf^d'x, (8) 

iOTT J 

where the variables are defined in Box 1, and where F^i, = A,y^^ — A^^^. By 
identifying qqq = T and gij = H6ij in a SSS field, Fio = Ei and Fij = eijkBk, 
one obtains Eq. (|^. Conversely, every theory within this class that satisfies 
Eq. d^) can have its electrodynamic equations cast into "metric" form. In a 
given non-metric theory, the functions T, H, e and /i will depend in general 
on the full gravitational environment, including the potential of the Earth, 
Sun and Galaxy, as well as on cosmological boundary conditions. Which of 
these factors has the most influence on a given experiment will depend on 
the nature of the experiment. 

Lightman and Lee then calculated explicitly the rate of fall of a "test" 
body made up of interacting charged particles, and found that the rate was 
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Box 1. The THeji Formalism 



1. Coordinate System and Conventions: 

= t = time coordinate associated with the static nature of the static 
spherically symmetric (sss) gravitational field; x = (x, y, z) = isotropic 
quasi-Cartesian spatial coordinates; spatial vector and gradient operations 
as in Cartesian space. 

2. Matter and Field Variables: 

• i^Oa = rest mass of particle a. 

• Ca = charge of particle a. 

• x'^{t) = world line of particle a. 

• = dx^^/dt = coordinate velocity of particle a. 

• = electromagnetic vector potential; E = VAQ — dA/dt , B = V x A 

3. Gravitational Potential: U{yi) 

4. Arbitrary Functions: 

T{U), H{U), e{U), ii{U); EEP is satisfied iff e = = {RlTfl'^ for ah U . 

5. Action: 

1 = -Y^moa [{T-Hvl)^/^dt + J2ea [ A^ixMdt 
+(87r)-i j {eE^ - fi-^B^)d^x . 

6. Non-Metric Parameters: 

To = -(^id/dU)ln[eiT/H)'%, 
Ao = -cUd/dU)ln[fi{T/H)y%, 
To = l-{TH-hfi)o, 

where cq = {Tq/Hq)^/^ and subscript "0" refers to a chosen point in space. 
If EEP is satisfied, Tq = Aq = Tq = 0. 
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independent of the internal electromagnetic structure of the body (WEP) 
if and only if Eq. was satisfied. In other words WEP EEP and 
Schiff 's conjecture was verified, at least within the restrictions built into the 
formalism. 

Certain combinations of the functions T, H, e and /i reflect different 
aspects of EEP. For instance, position or [/-dependence of either of the com- 
binations e(T/if)^/^ and /i(T/if)^/^ signals violations of LPI, the first combi- 
nation playing the role of the locally measured electric charge or fine structure 
constant. The "non-metric parameters" Fq and Aq (Box 1) are measures of 
such violations of EEP. Similarly, if the parameter Tq = 1 — {TH~^efi)o is 
non-zero anywhere, then violations of LLI will occur. This parameter is re- 
lated to the difference between the speed of light, c, and the limiting speed 
of material test particles, Cq, given by 

c=(eo/io)''/', Co = {To/H,f'\ (9) 

In many applications, by suitable definition of units, cq can be set equal to 
unity. If EEP is valid, Fq = Aq = Tq = everywhere. 

The rate of fall of a composite spherical test body of electromagnetically 
interacting particles then has the form 

a = (mp/m)V[/, (10) 

mp/m = 1 + (E|^/Mc2) [2Fo-^To] + (i?r/Mc^)[2Ao-^To] 

+ ..., (11) 

where E^^ and are the electrostatic and magnetostatic binding energies 
of the body, given by 



ab 



^ \ab 




where Vab = — ^b\, ^ab = (xa — ^b)/fab, and the angle brackets denote 
an expectation value of the enclosed operator for the system's internal state. 
Eotvos experiments place limits on the WEP- violating terms in Eq. ([Tl|) , 
and ultimately place limits on the non-metric parameters |Fo| < 2 x 10~^° 
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and |Ao| < 3 x 10^^. (We set Tq = because of very tight constraints on it 
from tests of LLI .) These hmits are sufficiently tight to rule out a number of 
non-metric theories of gravity thought previously to be viable (TEGP 2.6(f)). 

The Tife/i formalism also yields a gravitationally modified Dirac equation 
that can be used to determine the gravitational redshift experienced by a 
variety of atomic clocks. For the redshift parameter a [Eq. (^)], the results 
are (TEGP 2.6(c)): 

— SFq + Aq hydrogen hyperfine transition, H — Maser clock 
— |(3ro + Aq) electromagnetic mode in cavity, SCSO clock 



a 



— 2ro phonon mode in solid, principal transition in 

hydrogen. 

(14) 

The redshift is the standard one (a = 0), independently of the nature of 
the clock if and only if Fq = Aq = 0. Thus the Vessot-Levine rocket redshift 
experiment sets a limit on the parameter combination |3Fo — Ao| (Figure 
^); the null- redshift experiment comparing hydrogen-maser and SCSO clocks 
sets a limit on \au — olscso\ = flTo — Ao|. Alvarez and Mann ||^, ^, |, |, |^ 
extended the THen formalism to permit analysis of such effects as the Lamb 
shift, anomalous magnetic moments and non-baryonic effects, and placed 
interesting bounds on EEP violations. 



2.2.3 The Formalism 

The THefi formalism can also be applied to tests of local Lorentz invariance, 
but in this context it can be simplified. Since most such tests do not concern 
themselves with the spatial variation of the functions T, H, e, and /i, but 
rather with observations made in moving frames, we can treat them as spatial 
constants. Then by rescaling the time and space coordinates, the charges 
and the electromagnetic fields, we can put the action in Box 1 into the form 
(TEGP 2.6(a)). 

/ = -E^oa / {l-viy/'dt + Y^ea I A^«)<rft+(87r)-i / {E^-^B^)d'x, 

(15) 

where = Ho/ToeofiQ = (1 — Tq) ^. This amounts to using units in which 
the limiting speed Co of massive test particles is unity, and the speed of light 
is c. If c 7^ 1, LLI is violated; furthermore, the form of the action above 
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must be assumed to be valid only in some preferred universal rest frame. 
The natural candidate for such a frame is the rest frame of the microwave 
background. 

The electro dynamical equations which follow from Eq. ([ISD yield the 
behavior of rods and clocks, just as in the full THen formalism. For example, 
the length of a rod moving through the rest frame in a direction parallel to 
its length will be observed by a rest observer to be contracted relative to 
an identical rod perpendicular to the motion by a factor 1 — V"^ /2 + OiV^). 
Notice that c does not appear in this expression. The energy and momentum 
of an electromagnetically bound body which moves with velocity V relative 
to the rest frame are given by 

E = MR + ^MnV^ + ^6MyV'V\ 

= MrV' + 6Myv\ (16) 

where Mji = Mq — , Mq is the sum of the particle rest masses, is the 
electrostatic binding energy of the system ([Eq. (|l|)] with To^^HoeQ^ = 1), 
and 

SMy = -2{\- l)[tE^'S'^ + E^'''] , (17) 

where 

= -\ (E ^aCr-,' [^K,ni, - i^^^) ^ . (18) 
Note that (c~^ — 1) corresponds to the parameter 6 plotted in Figure ^. 



The electrodynamics given by Eq. ([TsD can also be quantized, so that 
we may treat the interaction of photons with atoms via perturbation theory. 
The energy of a photon is h times its frequency u, while its momentum is 
huj/c. Using this approach, one finds that the difference in round trip travel 
times of light along the two arms of the interferometer in the Michelson- 
Morley experiment is given by Lo(f^/c)(c~^ — 1). The experimental null 
result then leads to the bound on (c^^ — 1) shown on Figure 0. Similarly 
the anisotropy in energy levels is clearly illustrated by the tensorial terms 
in Eqs. (|l^) and (|l^); by evaluating E^^^-' for each nucleus in the various 
Hughes-Drever-type experiments and comparing with the experimental limits 
on energy differences, one obtains the extremely tight bounds also shown on 
Figure ||. 
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The behavior of moving atomic clocks can also be analysed in detail, and 
bounds on (c~^ — 1) can be placed using results from tests of time dilation and 
of the propagation of light. In some cases, it is advantageous to combine the 
(? framework with a "kinematical" viewpoint that treats a general class of 
boost transformations between moving frames. Such kinematical approaches 
have been discussed by Robertson, Mansouri and Sexl, and Will (see ||141|| ). 

For example, in the "JPL" experiment, in which the phases of two hy- 
drogen masers connected by a fiberoptic link were compared as a function 
of the Earth's orientation, the predicted phase difference as a function of 
direction is, to first order in V, the velocity of the Earth through the cosmic 
background, 

A0/0 ^ -^(1 - c2)(V ■ n - V ■ no) , (19) 

where = 27rz/L, v is the maser frequency, L = 21 km is the baseline, and 
where n and ng are unit vectors along the direction of propagation of the 
light, at a given time, and at the initial time of the experiment, respectively. 
The observed limit on a diurnal variation in the relative phase resulted in 
the bound |c~^ — 1| < 3 x 10~^. Tighter bounds were obtained from a "two- 
photon absorption" (TPA) experiment, and a 1960s series of "Mossbauer- 
rotor" experiments, which tested the isotropy of time dilation between a 
gamma ray emitter on the rim of a rotating disk and an absorber placed at 



the center ftL41 



2.3 EEP, Particle Physics, and the Search for New In- 
teractions 

In 1986, as a result of a detailed reanalysis of Eotvos' original data, Fischbach 
et al. |6^ suggested the existence of a fifth force of nature, with a strength 
of about a percent that of gravity, but with a range (as defined by the range 
A of a Yukawa potential, e^^^^ /r) of a few hundred meters. This proposal 
dovetailed with earlier hints of a deviation from the inverse-square law of 
Newtonian gravitation derived from measurements of the gravity profile down 
deep mines in Australia, and with ideas from particle physics suggesting the 
possible presence of very low-mass particles with gravitational-strength cou- 
plings. During the next four years numerous experiments looked for evidence 
of the fifth force by searching for composition-dependent differences in accel- 
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eration, with variants of the Eotvos experiment or with free-fall Galileo-type 
experiments. Although two early experiments reported positive evidence, the 
others all yielded null results. Over the range between one and lO'^ meters, 
the null experiments produced upper limits on the strength of a postulated 
fifth force between 10~^ and 10~^ of the strength of gravity. Interpreted as 
tests of WEP (corresponding to the limit of infinite-range forces), the results 
of two representative experiments from this period: the free-fall Galileo ex- 
periment, and the early Eot-Wash experiment, are shown in Figure |l[ At 
the same time, tests of the inverse-square law of gravity were carried out 
by comparing variations in gravity measurements up tall towers or down 
mines or boreholes with gravity variations predicted using the inverse square 
law together with Earth models and surface gravity data mathematically 
"continued" up the tower or down the hole. Despite early reports of anoma- 
lies, independent tower, borehole and seawater measurements now show no 
evidence of a deviation. Analyses of orbital data from planetary range mea- 
surements, lunar laser ranging, and laser tracking of the LAGEOS satellite 
verified the inverse-square law to parts in 10*^ over scales of 10^ to 10^ km, and 
to parts in 10^ over planetary scales of several astronomical units ||122||. The 



consensus at present is that there is no credible experimental evidence for a 
fifth force of nature. For reviews and bibhographies, see |^ ^ |], p.4CI|| . 

Nevertheless, theoretical evidence continues to mount that EEP is likely to 
be violated at some level, whether by quantum gravity effects, by effects aris- 
ing from string theory, or by hitherto undetected interactions, albeit at levels 
well below those that motivated the fifth-force searches. Roughly speaking, 
in addition to the pure Einsteinian gravitational interaction, which respects 
EEP, theories such as string theory predict other interactions which do not. 
In string theory, for example, the existence of such EEP- violating fields is as- 
sured, but the theory is not yet mature enough to enable calculation of their 
strength (relative to gravity), or their range (whether they are long range, 
like gravity, or short range, like the nuclear and weak interactions, and thus 
too short-range to be detectable). 

In one simple example, one can write the Lagrangian for the low-energy 
limit of string theory in the so-called "Einstein frame", in which the gravi- 
tational Lagrangian is purely general relativistic: 
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(20) 



where g^^ is the non-physical metric, R^y is the Ricci tensor derived from it, 
(/9 is a dilaton field, and G, V and M are functions of Lp. The Lagrangian 
includes that for the electromagnetic field F^^,, and that for particles, written 
in terms of Dirac spinors tl). This is not a metric representation because of 
the coupling of to matter via M(y9) and U{ip). A conformal transformation 
g^j,^ = F{(f)g^,^, 'ip = F{Lp)~^^'^ip, puts the Lagrangian in the form ("Jordan" 
frame) 
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One may choose F{ip) = const/M(y9)^ so that the particle Lagrangian takes 
the metric form (no explicit coupling to ip), but the electromagnetic La- 
grangian will still couple non-metrically to U{(f). The gravitational La- 
grangian here takes the form of a scalar-tensor theory (Sec. |3.3.2|) . But 
the non-metric electromagnetic term will, in general, produce violations of 



EEP. For examples of specific models, see [125, 49 



Thus, EEP and related tests are now viewed as ways to discover or place 
constraints on new physical interactions, or as a branch of "non-accelerator 
particle physics", searching for the possible imprints of high-energy parti- 
cle effects in the low-energy realm of gravity. Whether current or proposed 
experiments can actually probe these phenomena meaningfully is an open 
question at the moment, largely because of a dearth of firm theoretical pre- 
dictions. Despite this uncertainty, a number of experimental possibilities are 
being explored. 

Concepts for an equivalence principle experiment in space have been de- 
veloped. The project MICROSCOPE, designed to test WEP to lO'^^ has 
been approved by the French space agency CNES for a possible 2004 launch. 
Another, known as Satellite Test of the Equivalence Principle (STEP), is 
under consideration as a possible joint effort of NASA and the European 
Space Agency (ESA), with the goal of a 10"^*^ test. The gravitational red- 
shift could be improved to the 10~^° level using atomic clocks on board a 
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spacecraft which would travel to within four solar radii of the Sun. Labora- 
tory tests of the gravitational inverse square law at sub-millimeter scales are 
being developed as ways to search for new short-range interactions or for the 
existence of large extra dimensions; the challenge of these experiments is to 
distinguish gravitation-like interactions from electromagnetic and quantum 



mechanical (Casimir) effects p8 



3 Tests of Post-Newtonian Gravity 

3.1 Metric Theories of Gravity and the Strong Equiv- 
alence Principle 

3.1.1 Universal Coupling and the Metric Postulates 

The overwhelming empirical evidence supporting the Einstein equivalence 
principle, discussed in the previous section, supports the conclusion that the 
only theories of gravity that have a hope of being viable are metric theories, 
or possibly theories that are metric apart from possible weak or short-range 
non- metric couplings (as in string theory). Therefore for the remainder of this 
article, we shall turn our attention exclusively to metric theories of gravity, 
which assume that (i) there exists a symmetric metric, (ii) test bodies follow 
geodesies of the metric, and (iii) in local Lorentz frames, the non-gravitational 
laws of physics are those of special relativity. 

The property that all non-gravitational fields should couple in the same 
manner to a single gravitational field is sometimes called "universal cou- 
pling". Because of it, one can discuss the metric as a property of spacetime 
itself rather than as a field over spacetime. This is because its properties may 
be measured and studied using a variety of different experimental devices, 
composed of different non-gravitational fields and particles, and, because of 
universal coupling, the results will be independent of the device. Thus, for 
instance, the proper time between two events is a characteristic of spacetime 
and of the location of the events, not of the clocks used to measure it. 

Consequently, if EEP is valid, the non-gravitational laws of physics may 
be formulated by taking their special relativistic forms in terms of the Min- 
kowski metric rj and simply "going over" to new forms in terms of the curved 
spacetime metric g, using the mathematics of differential geometry. The 
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details of this "going over" can be found in standard textbooks (|^, |135|] ; 
TEGP 3.2) 



3.1.2 The Strong Equivalence Principle 

In any metric theory of gravity, matter and non-gravitational fields respond 
only to the spacetime metric g. In principle, however, there could exist other 
gravitational fields besides the metric, such as scalar fields, vector fields, and 
so on. If, by our strict definition of metric theory, matter does not couple 
to these fields, what can their role in gravitation theory be? Their role 
must be that of mediating the manner in which matter and non-gravitational 
fields generate gravitational fields and produce the metric; once determined, 
however, the metric alone acts back on the matter in the manner prescribed 
by EEP. 

What distinguishes one metric theory from another, therefore, is the num- 
ber and kind of gravitational fields it contains in addition to the metric, 
and the equations that determine the structure and evolution of these fields. 
From this viewpoint, one can divide all metric theories of gravity into two 
fundamental classes: "purely dynamical" and "prior-geometric". 

By "purely dynamical metric theory" we mean any metric theory whose 
gravitational fields have their structure and evolution determined by coupled 
partial differential field equations. In other words, the behavior of each field 
is influenced to some extent by a coupling to at least one of the other fields 
in the theory. By "prior geometric" theory, we mean any metric theory that 
contains "absolute elements" , fields or equations whose structure and evolu- 
tion are given a priori, and are independent of the structure and evolution 
of the other fields of the theory. These "absolute elements" typically include 
flat background metrics ry, cosmic time coordinates t, algebraic relationships 
among otherwise dynamical fields, such as g^j^i/ = h^j^i/ -\- k^ky, where h^^y and 
fc^ may be dynamical fields. 

General relativity is a purely dynamical theory since it contains only one 
gravitational field, the metric itself, and its structure and evolution are gov- 
erned by partial differential equations (Einstein's equations). Brans-Dicke 
theory and its generalizations are purely dynamical theories; the field equa- 
tion for the metric involves the scalar field (as well as the matter as source), 
and that for the scalar field involves the metric. Rosen's bimetric theory is 
a prior-geometric theory: it has a non-dynamical, Riemann-fiat background 
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metric, 77, and the field equations for the physical metric g involve rj. 

By discussing metric theories of gravity from this broad point of view, 
it is possible to draw some general conclusions about the nature of gravity 
in different metric theories, conclusions that arc reminiscent of the Einstein 
equivalence principle, but that are subsumed under the name "strong equiv- 
alence principle". 

Consider a local, freely falling frame in any metric theory of gravity. Let 
this frame be small enough that inhomogeneities in the external gravitational 
fields can be neglected throughout its volume. On the other hand, let the 
frame be large enough to encompass a system of gravitating matter and its 
associated gravitational fields. The system could be a star, a black hole, 
the solar system or a Cavendish experiment. Call this frame a "quasi- local 
Lorentz frame" . To determine the behavior of the system we must calculate 
the metric. The computation proceeds in two stages. First wc determine 
the external behavior of the metric and gravitational fields, thereby estab- 
lishing boundary values for the fields generated by the local system, at a 
boundary of the quasi-local frame "far" from the local system. Second, we 
solve for the fields generated by the local system. But because the metric 
is coupled directly or indirectly to the other fields of the theory, its struc- 
ture and evolution will be infiuenced by those fields, and in particular by 
the boundary values taken on by those fields far from the local system. This 
will be true even if we work in a coordinate system in which the asymp- 
totic form of Qui, in the boundary region between the local system and the 
external world is that of the Minkowski metric. Thus the gravitational en- 
vironment in which the local gravitating system resides can infiuence the 
metric generated by the local system via the boundary values of the auxil- 
iary fields. Consequently, the results of local gravitational experiments may 
depend on the location and velocity of the frame relative to the external 
environment. Of course, local non- gravitational experiments are unaffected 
since the gravitational fields they generate are assumed to be negligible, and 
since those experiments couple only to the metric, whose form can always be 
made locally Minkowskian at a given spacetime event. Local gravitational 
experiments might include Cavendish experiments, measurement of the ac- 
celeration of massive self-gravitating bodies, studies of the structure of stars 
and planets, or analyses of the periods of "gravitational clocks" . We can now 
make several statements about different kinds of metric theories. 

(i) A theory which contains only the metric g yields local gravitational 
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physics which is independent of the location and velocity of the local system. 
This follows from the fact that the only field coupling the local system to 
the environment is g, and it is always possible to find a coordinate system 
in which g takes the Minkowski form at the boundary between the local 
system and the external environment. Thus the asymptotic values of g^^ are 
constants independent of location, and are asymptotically Lorentz invariant, 
thus independent of velocity. General relativity is an example of such a 
theory. 

(ii) A theory which contains the metric g and dynamical scalar fields tpA 
yields local gravitational physics which may depend on the location of the 
frame but which is independent of the velocity of the frame. This follows 
from the asymptotic Lorentz invariance of the Minkowski metric and of the 
scalar fields, but now the asymptotic values of the scalar fields may depend 
on the location of the frame. An example is Brans-Dicke theory, where the 
asymptotic scalar field determines the effective value of the gravitational 
constant, which can thus vary as if varies. On the other hand, a form of 
velocity dependence in local physics can enter indirectly if the asymptotic 
values of the scalar field vary with time cosmologically. Then the rate of 
variation of the gravitational constant could depend on the velocity of the 
frame. 

(iii) A theory which contains the metric g and additional dynamical vector 
or tensor fields or prior-geometric fields yields local gravitational physics 
which may have both location and velocity-dependent effects. 

These ideas can be summarized in the strong equivalence principle (SEP), 
which states that 

1. WEP is vahd for self-gravitating bodies as well as for test bodies. 

2. The outcome of any local test experiment is independent of the velocity 
of the (freely faUing) apparatus. 

3. The outcome of any local test experiment is independent of where and 
when in the universe it is performed. 

The distinction between SEP and EEP is the inclusion of bodies with self- 
gravitational interactions (planets, stars) and of experiments involving grav- 
itational forces (Cavendish experiments, gravimeter measurements). Note 
that SEP contains EEP as the special case in which local gravitational forces 
are ignored. 
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The above discussion of the couphng of auxihary fields to local gravitat- 
ing systems indicates that if SEP is strictly valid, there must be one and only 
one gravitational field in the universe, the metric g. These arguments are 
only suggestive however, and no rigorous proof of this statement is available 
at present. Empirically it has been found that every metric theory other 
than GR introduces auxiliary gravitational fields, either dynamical or prior 
geometric, and thus predicts violations of SEP at some level (here we ig- 
nore quantum-theory inspired modifications to GR involving "i?^" terms). 
General relativity seems to be the only metric theory that embodies SEP 
completely. This lends some credence to the conjecture SEP General Rel- 
ativity. In Sec. we shall discuss experimental evidence for the validity 
of SEP. 



3.2 The Parametrized Post-Newtonian Formalism 

Despite the possible existence of long-range gravitational fields in addition to 
the metric in various metric theories of gravity, the postulates of those theo- 
ries demand that matter and non-gravitational fields be completely oblivious 
to them. The only gravitational field that enters the equations of motion 
is the metric g. The role of the other fields that a theory may contain can 
only be that of helping to generate the spacetime curvature associated with 
the metric. Matter may create these fields, and they plus the matter may 
generate the metric, but they cannot act back directly on the matter. Matter 
responds only to the metric. 

Thus the metric and the equations of motion for matter become the pri- 
mary entities for calculating observable effects, and all that distinguishes 
one metric theory from another is the particular way in which matter and 
possibly other gravitational fields generate the metric. 

The comparison of metric theories of gravity with each other and with 
experiment becomes particularly simple when one takes the slow-motion, 
weak-field limit. This approximation, known as the post-Newtonian limit, 
is sufficiently accurate to encompass most solar-system tests that can be 
performed in the foreseeable future. It turns out that, in this limit, the 
spacetime metric g predicted by nearly every metric theory of gravity has 
the same structure. It can be written as an expansion about the Minkowski 
metric {t]^^ = diag(— 1, 1, 1, 1)) in terms of dimensionless gravitational poten- 
tials of varying degrees of smallness. These potentials are constructed from 
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Tabic 2: The PPN Parameters and their significance (note that has been 
shown twice to indicate that it is a measure of two effects) 



the matter variables (Box 2) in imitation of the Newtonian gravitational 
potential 

C/(x, t) = j p(x', t) |x - ^\-^Sx' . (22) 

The "order of smallness" is determined according to the rules C/ ~ ~ 
H ~ p/ p ~ e, V* ~ \d/ dt\/\d/ dx\ ~ e-*^/^, and so on (we use units in which 
G = c = 1; see Box 2). 

A consistent post-Newtonian limit requires determination of 5^00 correct 
through O(e^), goi through 0(e^/^) and Qij through 0(e) (for details see TEGP 
4.1). The only way that one metric theory differs from another is in the nu- 
merical values of the coefficients that appear in front of the metric potentials. 
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The parametrized post-Newtonian (PPN ) formalism inserts parameters in 
place of these coefficients, parameters whose values depend on the theory 
under study. In the current version of the PPN formalism, summarized in 
Box 2, ten parameters are used, chosen in such a manner that they measure 
or indicate general properties of metric theories of gravity (Table |^). Under 
reasonable assumptions about the kinds of potentials that can be present at 
post-Newtonian order (basically only Poisson-like potentials), one finds that 
ten PPN parameters exhaust the possibilities. 

The parameters 7 and (3 are the usual Eddington-Robertson-Schiff param- 
eters used to describe the "classical" tests of GR, and are in some sense the 
most important; they are the only non-zero parameters in GR and scalar- 
tensor gravity. The parameter ^ is non-zero in any theory of gravity that 
predicts preferred- location effects such as a galaxy-induced anisotropy in the 
local gravitational constant Gl (also called "Whitehead" effects); ai, 02, 03 
measure whether or not the theory predicts post-Newtonian preferred-frame 
effects; as, Ci, C25 Cs; C4 measure whether or not the theory predicts viola- 
tions of global conservation laws for total momentum. Next to 7 and /3, the 
parameters ai and 02 occur most frequently with non-trivial null values. In 
Table ^ we show the values these parameters take (i) in GR, (ii) in any the- 
ory of gravity that possesses conservation laws for total momentum, called 
"semi-conservative" (any theory that is based on an invariant action principle 
is semi-conservative), and (iii) in any theory that in addition possesses six 
global conservation laws for angular momentum, called "fully conservative" 
(such theories automatically predict no post-Newtonian preferred-frame ef- 
fects). Semi-conservative theories have five free PPN parameters (7, /?, ^, 
ai, 02) while fully conservative theories have three (7, [3 , ^). 

The PPN formalism was pioneered by Kenneth Nordtvedt , who stud- 
ied the post-Newtonian metric of a system of gravitating point masses, ex- 
tending earlier work by Eddington, Robertson and Schiff (TEGP 4.2). A 
general and unified version of the PPN formalism was developed by Will and 
Nordtvedt. The canonical version, with conventions altered to be more in 
accord with standard textbooks such as [Q], is discussed in detail in TEGP, 
Chapter 4. Other versions of the PPN formalism have been developed to deal 
with point masses with charge, fiuid with anisotropic stresses, bodies with 
strong internal gravity, and post-post-Newtonian effects (TEGP 4.2, 14.2). 
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Box 2. The Pcirametrized Post-Newtonian Formalism 



1. Coordinate System: The framework uses a nearly globally Lorentz coor- 
dinate system in which the coordinates are {t, ^). Three-dimensional, 
Euclidean vector notation is used throughout. All coordinate arbitrariness 
("gauge freedom") has been removed by specialization of the coordinates to 
the standard PPN gauge (TEGP 4.2). Units are chosen so that G = c = 1, 
where G is the physically measured Newtonian constant far from the solar 
system. 

2. Matter Variables: 

• p = density of rest mass as measured in a local freely falling frame 

momentarily comoving with the gravitating matter. 

• = {dx^/dt) = coordinate velocity of the matter. 

• w'' = coordinate velocity of PPN coordinate system relative to the 
mean rest-frame of the universe. 

• p = pressure as measured in a local freely falling frame momentarily 
comoving with the matter. 

• n = internal energy per unit rest mass. It includes all forms of non- 
rcst-mass, non-gravitational energy, e.g. energy of compression and 
thermal energy. 

3. PPN Parameters: 

4. Metric: 

500 = -l + 2?7-2/3C/2-2e$;y + (27 + 2 + 03 + Ci-20$i 

+2(37 - 2/3 + 1 + C2 + 6^2 + 2(1 + C3)<^3 + 2(37 + 3C4 - 20$4 
-(Ci - 2^)^ - (ai - a2 - a3)w'^U - a2w''wWij + (2a3 - aijvf'Vi 

90i = -^(47 + 3 + ai-a2 + Ci-20^i-^(l + «2-Ci + 20W^i 

-^(ai - 2a2)w'U - a2W^Uij + 0(e^/2) 
Qij = {l + 2-fU + 0{e^))dij 
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Box 2. (continued) 



1. Metric Potentials: 



7 |x — x'l'^ V|x — x''| |x' — x"|/ 

J X — x' J X — x' 



$2 = l-P^d^x', $3=/^^dV, $4=/^_d3^' 

J |x — x'l J |x — x'l y |x — x'l 



P'[V'-(X-X')](X-X0^^3^, 



|x — x'P 

2. Stress-Energy Tensor (perfect fluid) 

r°° = p(l + n + 2^7) 

= p7;*(l + n + 7;2 + 2[/+p/p) 

= py'v^il + U + v"^ + 2U + p/p)+pS'^{l~2jU) 

3. Equations of Motion 

• Stressed Matter, T^".^ = 

• Test Bodies, (fx^'/dX^ + T^" ^:,{dx'' / d\){dx^ / dX) = 

• Maxwell's Equations, F^^'' = AnJ^" F^,, = A^./^ - Afj^.^ 
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3.3 Competing Theories of Gravity 

One of the important applications of the PPN formahsm is the comparison 
and classification of alternative metric theories of gravity. The population of 
viable theories has fluctuated over the years as new effects and tests have been 
discovered, largely through the use of the PPN framework, which eliminated 
many theories thought previously to be viable. The theory population has 
also fluctuated as new, potentially viable theories have been invented. 

In this article, we shall focus on general relativity and the general class of 
scalar-tensor modifications of it, of which the Jordan-Fierz-Brans-Dicke the- 
ory (Brans-Dicke, for short) is the classic example. The reasons are several- 
fold: 

• A full compendium of alternative theories is given in Chapter 5 of 
TEGP. 

• Many alternative metric theories developed during the 1970s and 1980s 
could be viewed as "straw-man" theories, invented to prove that such 
theories exist or to illustrate particular properties. Few of these could 
be regarded as well-motivated theories from the point of view, say, of 
field theory or particle physics. Examples are the vector-tensor theories 
studied by Will, Nordtvedt and Hellings. 

• A number of theories fall into the class of "prior-geometric" theories, 
with absolute elements such as a flat background metric in addition to 
the physical metric. Most of these theories predict "preferred-frame" 
effects, that have been tightly constrained by observations (see Sec. 
^.6.2[ ). An example is Rosen's bimetric theory. 



A large number of alternative theories of gravity predict gravitational- 
wave emission substantially different from that of general relativity, in 
strong disagreement with observations of the binary pulsar (see Sec. 

Scalar-tensor modifications of GR have recently become very popular 
in unification schemes such as string theory, and in cosmological model 
building. Because the scalar fields are generally massive, the potentials 
in the post-Newtonian limit will be modified by Yukawa-like terms. 
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3.3.1 General Relativity 



The metric g is the sole dynamical field and the theory contains no arbitrary 
functions or parameters, apart from the value of the Newtonian coupling 
constant G, which is measurable in laboratory experiments. Throughout this 
article, we ignore the cosmological constant A. Although A has significance 
for quantum field theory, quantum gravity, and cosmology, on the scale of 
the solar-system or of stellar systems, its effects are negligible, for values of 
A corresponding to a cosmological closure density. 

The field equations of GR are derivable from an invariant action principle 
61 = 0, where 

/ = (167rG)-i J R{-gf'H^x + g^,) , (23) 

where R is the Ricci scalar, and is the matter action, which depends on 
matter fields ipm universally coupled to the metric g. By varying the action 
with respect to g^^, we obtain the field equations 

G^u = Rfj,u — -g^J.uR = SnCTfj^i, , (24) 

where T^i, is the matter energy-momentum tensor. General covariance of the 
matter action implies the equations of motion T^^^-u = 0; varying with 
respect to ipM yields the matter field equations. By virtue of the absence of 
prior-geometric elements, the equations of motion are also a consequence of 
the field equations via the Bianchi identities C^^.^y = 0. 

The general procedure for deriving the post-Newtonian limit is spelled 
out in TEGP 5.1, and is described in detail for GR in TEGP 5.2. The PPN 
parameter values are listed in Table ^ 



3.3.2 Scalar- Tensor Theories 



These theories contain the metric g, a scalar field (p, a potential function 
V{(p), and a coupling function A{(f) (generalizations to more than one scalar 
field have also been carried out PT|). For some purposes, the action is conve- 
niently written in a non-metric representation, sometimes denoted the "Ein- 
stein frame" , in which the gravitational action looks exactly like that of GR: 



/= (167rG)-^ / [R 



(25) 
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Table 3: Metric Theories and Their PPN Parameter Values {a^ = Q = for 
all cases) 



where R = g^^R^v is the Ricci scalar of the "Einstein" metric g^y. (Apart 
from the scalar potential term V{ip)^ this corresponds to Eq. ( PO] ) with 
G{ip) = {4:nG)~^, U{ip) = 1, and M{(p) oc ^(v^).) This representation is a 
"non-metric" one because the matter fields ipm couple to a combination of ip 
and gfj^i,. Despite appearances, however, it is a metric theory, because it can 
be put into a metric representation by identifying the "physical metric" 

g^y = A^{ip)g^y . (26) 

The action can then be rewritten in the metric form 

/ = (167rG)-i J [<PR - (j)-M<P)9^''d^<Pd.<p - (P'V]{-gy/'d^x + lUi^^n, 9,u) , 

(27) 

where 

= A{^)-\ 
3 + 2lj{4>) = a(v9)-^ 

a{ip) = dilnA{ip))/dip. (28) 

The Einstein frame is useful for discussing general characteristics of such 
theories, and for some cosmological applications, while the metric represen- 
tation is most useful for calculating observable effects. The field equations. 
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post-Newtonian limit and PPN parameters are discussed in TEGP 5.3, and 
the values of the PPN parameters are listed in Table ^ 
The parameters that enter the post-Newtonian limit are 

u = cj(0o) A = [(c/cu/rf0)(3 + 2cj)-2(4 + 2uj)-\, , (29) 

where 0o is the value of today far from the system being studied, as de- 
termined by appropriate cosmological boundary conditions. The following 
formula is also useful: 1/(2 + uj) = 20^/ {1 + a^). In Brans-Dicke theory 
(cij((/)) = constant), the larger the value of cj, the smaller the effects of the 
scalar field, and in the limit oj ^ oo (ao — > 0), the theory becomes indis- 
tinguishable from GR in all its predictions. In more general theories, the 
function cij(0) could have the property that, at the present epoch, and in 
weak-field situations, the value of the scalar field 0o is such that, uo is very 
large and A is very small (theory almost identical to GR today), but that 
for past or future values of 0, or in strong-field regions such as the interiors 
of neutron stars, uj and A could take on values that would lead to significant 
differences from GR. Indeed, Damour and Nordtvedt have shown that in 
such general scalar-tensor theories, GR is a natural "attractor" : regardless 
of how different the theory may be from GR in the early universe (apart 
from special cases), cosmological evolution naturally drives the fields toward 
small values of the function a, thence to large uj. Estimates of the expected 
relic deviations from GR today in such theories depend on the cosmological 
model, but range from 10~^ to a few times 10~^ for 1— 7 (|^, 0). 

Scalar fields coupled to gravity or matter are also ubiquitous in particle- 
physics-inspired models of unification, such as string theory. In some models, 
the coupling to matter may lead to violations of WEP, which are tested 
by Eotvos-type experiments. In many models the scalar field is massive; 
if the Compton wavelength is of macroscopic scale, its effects are those of 
a "fifth force". Only if the theory can be cast as a metric theory with 
a scalar field of infinite range or of range long compared to the scale of 
the system in question (solar system) can the PPN framework be strictly 
applied. If the mass of the scalar field is sufficiently large that its range is 
microscopic, then, on solar-system scales, the scalar field is suppressed, and 
the theory is essentially equivalent to general relativity. This is the case, for 
example in the "oscillating-G" models of Accetta, Steinhardt and Will (see 
[ |12CI|| ), in which the potential function V{ip) contains both quadratic (mass) 
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and quartic (self-interaction) terms, causing the scalar field to oscillate (the 
initial amplitude of oscillation is provided by an inflationary epoch); high- 
frequency oscillations in the "effective" Newtonian constant Gcff = G/cf) = 
GA{^f then result. The energy density in the oscillating scalar field can be 
enough to provide a cosmological closure density without resorting to dark 
matter, yet the value of uj today is so large that the theory's local predictions 
are experimentally indistinguishable from GR. In other models, explored by 
Damour and Esposito-Farese [|^, non-linear scalar-field couplings can lead 



to "spontaneous scalarization" inside strong-field objects such as neutron 
stars, leading to large deviations from GR, even in the limit of very large uo. 



3.4 Tests of the Parameter 7 

With the PPN formalism in hand, we are now ready to confront gravitation 
theories with the results of solar-system experiments. In this section we focus 
on tests of the parameter 7, consisting of the deflection of light and the time 
delay of light. 



3.4.1 The Deflection of Light 

A light ray (or photon) which passes the Sun at a distance d is deflected by 
an angle 

59 = ^(1 + 7)(4m0/c?)[(l + cos<l>)/2] (30) 

(TEGP 7.1), where itlq is the mass of the Sun and $ is the angle between 
the Earth-Sun line and the incoming direction of the photon (Figure ^). For 
a grazing ray, d ^ dQ, $ ^ 0, and 

56 ^ i(l + 7)l/'75, (31) 

independent of the frequency of light. Another, more useful expression gives 
the change in the relative angular separation between an observed source of 
light and a nearby reference source as both rays pass near the Sun: 



^^=^(1 + 7) 



4mQ Atuq / 1 + cos 

— 1- cos X + —r- 
d dr 



(32) 



where d and dr are the distances of closest approach of the source and ref- 
erence rays respectively, $r is the angular separation between the Sun and 
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Reference Source ft^ 




Figure 4: Geometry of light deflection measurements. 

the reference source, and x is the angle between the Sun-source and the Sun- 
reference directions, projected on the plane of the sky (Figure ^). Thus, for 
example, the relative angular separation between the two sources may vary 
if the line of sight of one of them passes near the Sun [d ~ i?0, dr ^ d, x 
varying with time). 

It is interesting to note that the classic derivations of the deflection of light 
that use only the principle of equivalence or the corpuscular theory of light 
yield only the "1/2" part of the coefficient in front of the expression in Eq. 
(|30D. But the result of these calculations is the deflection of light relative to 
local straight lines, as established for example by rigid rods; however, because 
of space curvature around the Sun, determined by the PPN parameter 7, 
local straight lines are bent relative to asymptotic straight lines far from the 
Sun by just enough to yield the remaining factor "7/2". The first factor 
"1/2" holds in any metric theory, the second "7/2" varies from theory to 
theory. Thus, calculations that purport to derive the full deflection using the 
equivalence principle alone are incorrect. 

The prediction of the full bending of light by the Sun was one of the 
great successes of Einstein's GR. Eddington's conflrmation of the bending 
of optical starlight observed during a solar eclipse in the flrst days following 
World War I helped make Einstein famous. However, the experiments of 
Eddington and his co-workers had only 30 percent accuracy, and succeeding 
experiments were not much better: the results were scattered between one 
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half and twice the Einstein value (Figure and the accuracies were low. 

However, the development of VLSI, very-long-baseline radio interferom- 
etry, produced greatly improved determinations of the deflection of light. 
These techniques now have the capability of measuring angular separations 
and changes in angles as small as 100 microarcseconds. Early measurements 
took advantage of a series of heavenly coincidences: each year, groups of 
strong quasistellar radio sources pass very close to the Sun (as seen from 
the Earth), including the group 3C273, 3C279, and 3C48, and the group 
0111+02, 0119+11 and 0116+08. As the Earth moves in its orbit, changing 
the lines of sight of the quasars relative to the Sun, the angular separation 59 
between pairs of quasars varies [Eq. (p2D ]. The time variation in the quanti- 
ties dr, X ci-nd $r in Eq- ( |52D is determined using an accurate ephemeris for 
the Earth and initial directions for the quasars, and the resulting prediction 
for 59 as a function of time is used as a basis for a least-squares fit of the 
measured 59, with one of the fitted parameters being the coefficient |(1 + 7). 
A number of measurements of this kind over the period 1969-1975 yielded 
an accurate determination of the coefficient ^(1 + 7) which has the value 
unity in GR. A 1995 VLB! measurement using 3C273 and 3C279 yielded 
(1 + 7)/2 = 0.9996 ± 0.0017 [0. 

A recent series of transcontinental and intercontinental VLBI quasar and 
radio galaxy observations made primarily to monitor the Earth's rotation 
( "VLBI " in Figure ^) was sensitive to the deflection of light over almost 
the entire celestial sphere (at 90° from the Sun, the deflection is still 4 milli- 
arcseconds). A recent analysis of over 2 million VLBI observations yielded 
(1 + 7)/2 = 0.99992 ± 0.00014 p[. Analysis of observations made by the 
Hipparcos optical astrometry satellite yielded a test at the level of 0.3 percent 
[|65|] . A VLBI measurement of the deflection of light by Jupiter was reported; 
the predicted deflection of about 300 microarcseconds was seen with about 
50 percent accuracy [129]. The results of light-deflection measurements are 
summarized in Figure ^. 



3.4.2 The Time Delay of Light 

A radar signal sent across the solar system past the Sun to a planet or satellite 
and returned to the Earth suffers an additional non-Newtonian delay in its 
round-trip travel time, given by (see Figure ^ 

5t = 2(1 + 7)m0 ln[(re + ■ n) (r^ - Xe ■ rv)/d^] (33) 
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Figure 5: Measurements of the coefficient (1 + 7)/2 from light deflection and 
time delay measurements. General relativity value is unity. Arrows denote 
anomalously large values from early eclipse expeditions. Shapiro time-delay 
measurements using Viking spacecraft yielded agreement with GR to 0.1 
percent, and VLBI light deflection measurements have reached 0.02 percent. 
Hipparcos denotes the optical astrometry satellite, which has reached 0.1 
percent. 



(TEGP 7.2). For a ray which passes close to the Sun, 

i(l + 7)[240 - 201n(rfVr)]/is, (34) 

where d is the distance of closest approach of the ray in solar radii, and r is 
the distance of the planet or satellite from the Sun, in astronomical units. 

In the two decades following Irwin Shapiro's 1964 discovery of this ef- 
fect as a theoretical consequence of general relativity, several high-precision 
measurements were made using radar ranging to targets passing through su- 
perior conjunction. Since one does not have access to a "Newtonian" signal 
against which to compare the round-trip travel time of the observed signal, 
it is necessary to do a differential measurement of the variations in round- 
trip travel times as the target passes through superior conjunction, and to 
look for the logarithmic behavior of Eq. (|3^ ) . In order to do this accurately 
however, one must take into account the variations in round-trip travel time 
due to the orbital motion of the target relative to the Earth. This is done by 
using radar-ranging (and possibly other) data on the target taken when it is 
far from superior conjunction (z.e., when the time-delay term is negligible) 
to determine an accurate ephemeris for the target, using the ephemeris to 
predict the PPN coordinate trajectory Xe(t) near superior conjunction, then 
combining that trajectory with the trajectory of the Earth X0(t) to deter- 
mine the Newtonian round-trip time and the logarithmic term in Eq. (|3^). 
The resulting predicted round-trip travel times in terms of the unknown co- 
efficient |(1 + 7) are then fit to the measured travel times using the method 
of least-squares, and an estimate obtained for ^(1 + 7). 

The targets employed included planets, such as Mercury or Venus, used 
as a passive refiectors of the radar signals ("passive radar"); and artificial 
satellites, such as Mariners 6 and 7, Voyager 2, and the Viking Mars lan- 
ders and orbiters, used as active retransmitters of the radar signals ( "active 
radar'). 

The results for the coefficient |(1 -f- 7) of all radar time-delay measure- 
ments performed to date (including a measurement of the one-way time delay 
of signals from the millisecond pulsar PSR 1937-1-21) are shown in Figure ^ 
(see TEGP 7.2 for discussion and references). The Viking experiment re- 
sulted in a 0.1 percent measurement ||111|| . 

From the results of VLBI light-defiection experiments, we can conclude 
that the coefficient |(1 -|- 7) must be within at most 0.014 percent of unity. 
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Parameter Effect 


Limit 


Remarks 


7-1 


time delay 


2 X 10-3 


Viking ranging 




light denection 


3 X 10-^ 


VLBI 


0-1 


perihelion shift 


3 X 10-3 


J2 = 10-^ from helioseismology 




Nordtvedt effect 


6 X 10-^ 


77 = 4/3 — 7 — 3 assumed 




Earth tides 


10-3 


gravimeter data 


ai 


orbital polarization 


10 ^ 


Lunar laser ranging 






2 X 10-^ 


PSR J2317+1439 




spin precession 


4 X 10-^ 


solar alignment with ecliptic 




pulsar acceleration 


2 X 10-20 


pulsar P statistics 




Nordtvedt effect 


10-3 


lunar laser ranging 


Ci 




2 X 10-2 


combined PPN bounds 


C2 


binary acceleration 


4 X 10-^ 


Pp for PSR 1913+16 


Cs 


Newton's 3rd law 


io-« 


Lunar acceleration 


Ci 






not in(lo])on(lont 


^Here 77 = 


4/3-7-3-10^/3- 


QJi — 2q;2/3 


- 2C1/3 - C2/3 



Table 4: Current Limits on the PPN Parameters 



Scalar-tensor theories must have u > 3500 to be compatible with this con- 
straint. 

3.5 The Perihelion Shift of Mercury 

The explanation of the anomalous perihelion shift of Mercury's orbit was an- 
other of the triumphs of GR. This had been an unsolved problem in celestial 
mechanics for over half a century, since the announcement by Le Verrier in 
1859 that, after the perturbing effects of the planets on Mercury's orbit had 
been accounted for, and after the effect of the precession of the equinoxes 
on the astronomical coordinate system had been subtracted, there remained 
in the data an unexplained advance in the perihelion of Mercury. The mod- 
ern value for this discrepancy is 43 arcseconds per century. A number of 
ad hoc proposals were made in an attempt to account for this excess, in- 
cluding, among others, the existence of new planet Vulcan near the Sun, 
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a ring of planetoids, a solar quadrupole moment and a deviation from the 
inverse-square law of gravitation, but none was successful. General relativity 
accounted for the anomalous shift in a natural way without disturbing the 
agreement with other planetary observations. 

The predicted advance, AcD , per orbit, including both relativistic PPN 
contributions and the Newtonian contribution resulting from a possible solar 
quadrupole moment, is given by 

Au = {6nm/p) [1 (2 + 27 - /5) + ^ (2ai - ^2 + "3 + 2(2)/^/^ + J2 {R^/2mp)] , 

6 D 

(35) 

where m = mi + m2 and yU = mim2/m are the total mass and reduced mass 
of the two-body system respectively; p = a(l — e^) is the semi-latus rectum 
of the orbit, with a the semi-major axis and e the eccentricity; R is the mean 
radius of the oblate body; and J2 is a dimensionless measure of its quadrupole 
moment, given by J2 = {C — A)/miR'^, where C and A are the moments 
of inertia about the body's rotation and equatorial axes, respectively (for 
details of the derivation see TEGP 7.3). We have ignored preferred-frame 
and galaxy-induced contributions to Au; these are discussed in TEGP 8.3. 

The first term in Eq. (P3| ) is the classical relativistic perihelion shift, which 
depends upon the PPN parameters 7 and /?. The second term depends upon 
the ratio of the masses of the two bodies; it is zero in any fully conservative 
theory of gravity (ai = 02 = 03 = C2 = 0); it is also negligible for Mercury, 
since /x/m ~ mMerc/^& ~ 2 x 10"''. We shall drop this term henceforth. 
The third term depends upon the solar quadrupole moment J2- For a Sun 
that rotates uniformly with its observed surface angular velocity, so that the 
quadrupole moment is produced by centrifugal flattening, one may estimate 
J2 to be ~ 1 X 10~^. This actually agrees reasonably well with values inferred 
from rotating solar models that are in accord with observations of the normal 
modes of solar oscillations (helioseismology) . Substituting standard orbital 
elements and physical constants for Mercury and the Sun we obtain the rate 
of perihelion shift u, in seconds of arc per century. 



u = 42."98 



^(2 + 27-/5) + 3 X 10-^(J2/10-') 



(36) 



Now, the measured perihelion shift of Mercury is known accurately: after the 
perturbing effects of the other planets have been accounted for, the excess 
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shift is known to about 0.1 percent from radar observations of Mercury be- 
tween 1966 and 1990 [ 116|| . Analysis of data taken since 1990 could improve 
the accuracy. The solar oblateness effect is smaller than the observational 
error, so we obtain the PPN bound \2'y — (3 — 1\ <3x 10~^. 



3.6 Tests of the Strong Equivalence Principle 

The next class of solar-system experiments that test relativistic gravitational 
effects may be called tests of the strong equivalence principle (SEP). In Sec. 
3.1.2, we pointed out that many metric theories of gravity (perhaps all ex- 
cept GR) can be expected to violate one or more aspects of SEP. Among 
the testable violations of SEP are a violation of the weak equivalence prin- 
ciple for gravitating bodies that leads to perturbations in the Earth-Moon 
orbit; preferred-location and preferred-frame effects in the locally measured 
gravitational constant that could produce observable geophysical effects; and 
possible variations in the gravitational constant over cosmological timescales. 



3.6.1 The Nordtvedt Effect and the Lunar Eotvos Experiment 

In a pioneering calculation using his early form of the PPN formalism, Nord- 



tvedt I93] showed that many metric theories of gravity predict that massive 



bodies violate the weak equivalence principle - that is, fall with different 



accelerations depending on their gravitational self-energy. Dicke |Q linked 
such an effect to the possibility of a spatially varying gravitational constant, 
in theories such as scalar-tensor gravity. For a spherically symmetric body, 
the acceleration from rest in an external gravitational potential U has the 
form 

a = {mp/m)VU , 
mp/m = l — ri[Eg/m), 

10 2 2 1 

= 4/3-7-3-ye-«i+ 3«2--Ci- 3C2, (37) 

where Eg is the negative of the gravitational self-energy of the body {Eg > 0). 
This violation of the massive-body equivalence principle is known as the 
"Nordtvedt effect" . The effect is absent in GR (77 = 0) but present in scalar- 
tensor theory {rj = l/(2+co')+4A). The existence of the Nordtvedt effect does 
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not violate the results of laboratory Eotvos experiments, since for laboratory- 
sized objects, Eg/m < 10"^'^, far below the sensitivity of current or future 
experiments. However, for astronomical bodies, Eg/m may be significant 
(10-^ for the Sun, 10"^ for Jupiter, 4.6 x 10"^° for the Earth, 0.2 x 10-^° 
for the Moon). If the Nordtvedt effect is present (^7 7^ 0) then the Earth 
should fall toward the Sun with a slightly different acceleration than the 
Moon. This perturbation in the Earth-Moon orbit leads to a polarization of 
the orbit that is directed toward the Sun as it moves around the Earth- Moon 
system, as seen from Earth. This polarization represents a perturbation in 
the Earth-Moon distance of the form 



5r = l?>.lrjco?,{ujQ — uJs)t m, (38) 

where o^o and Us are the angular frequencies of the orbits of the Moon and 
Sun around the Earth (see TEGP 8.1, for detailed derivations and references; 
for improved calculations of the numerical coefficient, see ||102| , ^). 



Since August 1969, when the first successful acquisition was made of a 
laser signal reflected from the Apollo 11 retrorefiector on the Moon, the lu- 
nar laser-ranging experiment (LURE) has made regular measurements of the 
round-trip travel times of laser pulses between a network of observatories 
and the lunar retrorefiectors, with accuracies that are approaching 50 ps (1 
cm) . These measurements are fit using the method of least-squares to a the- 
oretical model for the lunar motion that takes into account perturbations 
due to the Sun and the other planets, tidal interactions, and post-Newtonian 
gravitational effects. The predicted round-trip travel times between retrore- 
fiector and telescope also take into account the librations of the Moon, the 
orientation of the Earth, the location of the observatory, and atmospheric 
effects on the signal propagation. The "Nordtvedt" parameter r] along with 
several other important parameters of the model are then estimated in the 
least-squares method. 

Several independent analyses of the data found no evidence, within ex- 
perimental uncertainty, for the Nordtvedt effect (for recent results see 
|153| , p6|). Their results can be summarized by the bound \ri\ < 0.001. These 
results represent a limit on a possible violation of WEP for massive bodies 
of 5 parts in 10^^ (compare Figure]^). For Brans-Dicke theory, these results 
force a lower hmit on the coupling constant uj of 1000. Note that, at this 
level of precision, one cannot regard the results of lunar laser ranging as a 
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"clean" test of SEP until one eliminates the possibility of a compensating vi- 
olation of WEP for the two bodies, because the chemical compositions of the 
Earth and Moon differ in the relative fractions of iron and silicates. To this 
end, the Eot-Wash group carried out an improved test of WEP for labora- 
tory bodies whose chemical compositions mimic that of the Earth and Moon. 
The resulting bound of four parts in 10^'^ reduces the ambiguity in the 
Lunar laser ranging bound, and establishes the firm limit on the universality 
of acceleration of gravitational binding energy at the level of 77 < 1.3 x 10"'^. 

In GR, the Nordtvedt effect vanishes; at the level of several centimeters 
and below, a number of non-null general relativistic effects should be present 



3.6.2 Preferred-Frame and Preferred-Location Effects 

Some theories of gravity violate SEP by predicting that the outcomes of local 
gravitational experiments may depend on the velocity of the laboratory rela- 
tive to the mean rest frame of the universe (preferred-frame effects) or on the 
location of the laboratory relative to a nearby gravitating body (preferred- 
location effects). In the post-Newtonian limit, preferred-frame effects are 
governed by the values of the PPN parameters ai, 02, and 03, and some 
preferred-location effects are governed by ^ (see Table |^). 

The most important such effects are variations and anisotropics in the 
locally-measured value of the gravitational constant, which lead to anoma- 
lous Earth tides and variations in the Earth's rotation rate; anomalous con- 
tributions to the orbital dynamics of planets and the Moon; self-accelerations 
of pulsars, and anomalous torques on the Sun that would cause its spin axis 
to be randomly oriented relative to the ecliptic (see TEGP 8.2, 8.3, 9.3 and 
14.3(c)). An improved bound on 03 of 2 x 10~^° from the period derivatives 
of 20 millisecond pulsars was reported in [|11|; improved bounds on ai were 



achieved using lunar laser ranging data |^, and using observations of the 
circular binary orbit of the pulsar J2317+1439 (JT^j). Negative searches for 
these effects have produced strong constraints on the PPN parameters (Table 
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Method 



6-/^(10-12 yr-i) 



Lunar Laser Ranging ± 8 

3±5 

Viking Radar 2 ± 4 

-2± 10 

Binary Pulsar^ 11 ± 11 

Pulsar PSR 0655+64^ < 55 

1 Bounds dependent upon theory of gravity in 
strong-field regime and on neutron star equation 
of state. 

Table 5: Constancy of the Gravitational Constant 



3.6.3 Constancy of the Newtonian Gravitational Constant 

Most theories of gravity that violate SEP predict that the locally measured 
Newtonian gravitational constant may vary with time as the universe evolves. 
For the scalar-tensor theories listed in Table the predictions for G/G can 
be written in terms of time derivatives of the asymptotic scalar field. Where 
G does change with cosmic evolution, its rate of variation should be of the 
order of the expansion rate of the universe, i.e., G/G ~ Hq, where Hq is the 
Hubble expansion parameter and is given by Hq = lOO/i km Mpc"^ = 
h X 10^1'^ yr^i, where current observations of the expansion of the universe 
give h ^ 0.7. 

Several observational constraints can be placed on G/G using methods 
that include studies of the evolution of the Sun, observations of lunar occul- 
tations (including analyses of ancient eclipse data), lunar laser-ranging mea- 
surements, planetary radar-ranging measurements, and pulsar timing data. 
Laboratory experiments may one day lead to interesting limits (for review 
and references to past work see TEGP 8.4 and 14.3(c)). Recent results are 
shown in Table ^ 

The best limits on G/G still come from ranging measurements to the 



Viking landers and Lunar laser ranging measurements [^, |1 53| , [96[] . It has 



been suggested that radar observations of a Mercury orbiter over a two-year 
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mission (30 cm accuracy in range) could yield A{G/G) ~ 10~^^ yr"^- 

Although bounds on G/G from solar-system measurements can be cor- 
rectly obtained in a phenomenological manner through the simple expedient 
of replacing G hj Gq + Go(t — tg) in Newton's equations of motion, the same 
does not hold true for pulsar and binary pulsar timing measurements. The 
reason is that, in theories of gravity that violate SEP, such as scalar-tensor 
theories, the "mass" and moment of inertia of a gravitationally bound body 
may vary with variation in G. Because neutron stars are highly relativistic, 
the fractional variation in these quantities can be comparable to AG/G, the 
precise variation depending both on the equation of state of neutron star 
matter and on the theory of gravity in the strong-field regime. The variation 
in the moment of inertia affects the spin rate of the pulsar, while the varia- 
tion in the mass can affect the orbital period in a manner that can subtract 
from the direct effect of a variation in G, given by Pt/Pb = —2G/G | |1(J1| ]. 
Thus, the bounds quoted in Table |^ for the binary pulsar PSR 1913-1-16 |^ 



and the pulsar PSR 0655+64 ||69l are theory-dependent and must be treated 
as merely suggestive. 

3.7 Other Tests of Post-Newtonian Gravity 

3.7.1 Tests of Post-Newtonian Conservation Laws 

Of the five "conservation law" PPN parameters Ci, C2, Cs; C4; cind ^3, only 
three, (2, Ca ^ind a^, have been constrained directly with any precision; ^1 is 
constrained indirectly through its appearance in the Nordtvedt effect param- 
eter 1], Eq. (^). There is strong theoretical evidence that C4, which is related 
to the gravity generated by fluid pressure, is not really an independent pa- 
rameter - in any reasonable theory of gravity there should be a connection 
between the gravity produced by kinetic energy (pf^), internal energy (pll), 



and pressure (p). From such considerations, there follows [^3^ the additional 
theoretical constraint 

6C4 = 3«3 + 2Ci-3C3. (39) 

A non-zero value for any of these parameters would result in a violation of 
conservation of momentum, or of Newton's third law in gravitating systems. 
An alternative statement of Newton's third law for gravitating systems is 
that the "active gravitational mass", that is the mass that determines the 
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gravitational potential exhibited by a body, should equal the "passive gravita- 
tional mass" , the mass that determines the force on a body in a gravitational 
field. Such an equality guarantees the equality of action and reaction and of 
conservation of momentum, at least in the Newtonian limit. 

A classic test of Newton's third law for gravitating systems was carried 
out in 1968 by Kreuzer, in which the gravitational attraction of fiuorine and 
bromine were compared to a precision of 5 parts in 10^. 

A remarkable planetary test was reported by Bartlett and van Buren 
[ |T0[| . They noted that current understanding of the structure of the Moon 



involves an iron-rich, aluminum-poor mantle whose center of mass is offset 
about 10 km from the center of mass of an aluminum-rich, iron-poor crust. 
The direction of offset is toward the Earth, about 14° to the east of the 
Earth-Moon line. Such a model accounts for the basaltic maria which face 
the Earth, and the aluminum-rich highlands on the Moon's far side, and for 
a 2 km offset between the observed center of mass and center of figure for 
the Moon. Because of this asymmetry, a violation of Newton's third law for 
aluminum and iron would result in a momentum non-conserving self-force on 
the Moon, whose component along the orbital direction would contribute to 
the secular acceleration of the lunar orbit. Improved knowledge of the lunar 
orbit through lunar laser ranging, and a better understanding of tidal effects 
in the Earth- Moon system (which also contribute to the secular acceleration) 
through satellite data, severely limit any anomalous secular acceleration, with 
the resulting limit 



{mA/mp)A\ - {mA/mp)Ye 



{niA/mp) 



Fe 



< 4 X 10"^^ . (40) 



According to the PPN formalism, in a theory of gravity that violates conser- 
vation of momentum, but that obeys the constraint of Eq. (^91), the electro- 
static binding energy E^, of an atomic nucleus could make a contribution to 
the ratio of active to passive mass of the form 

mA = mp + ^(sEe/c^ ■ (41) 

The resulting limit on from the lunar experiment is < 1 x 10^*^ (TEGP 
9.2, 14.3(d)). 
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Another consequence of a violation of conservation of momentum is a 
self-acceleration of the center of mass of a binary stellar system, given by 



1,, miiSm e 
acM = - C2 + "3 — jz iT^np , 42) 

where 5m = mi — m2, a is the semi-major axis, and np is a unit vector 
directed from the center of mass to the point of periastron of mi (TEGP 
9.3). A consequence of this acceleration would be non- vanishing values for 
d?P/dt^, where P denotes the period of any intrinsic process in the system 
(orbit, spectra, pulsar periods). The observed upper limit on (PPp/dt^ of 
the binary pulsar PSR 1913-1-16 places a strong constraint on such an effect, 
resulting in the bound \a?j -|- C2I < 4 x 10~^. Since has already been 
constrained to be much less than this (Table §), we obtain a strong bound 
on C,2 alone |142|. 



3.7.2 Geodetic Precession 

A gyroscope moving through curved spacetime suffers a precession of its axis 
given by 

dS/dr = f2GxS, nG=(7 + ^)vx Vf/, (43) 

where v is the velocity of the gyroscope, and U is the Newtonian gravita- 
tional potential of the source (TEGP 9.1). The Earth- Moon system can be 
considered as a "gyroscope" , with its axis perpendicular to the orbital plane. 
The predicted precession is about 2 arcseconds per century, an effect first 
calculated by de Sitter. This effect has been measured to about 0.7 percent 



using Lunar laser ranging data |55, 153 



For a gyroscope orbiting the Earth, the precession is about 8 arcseconds 
per year. The Stanford Gyroscope Experiment has as one of its goals the 
measurement of this effect to 5 x 10^^ (see below); if achieved, this would 
substantially improve the accuracy of the parameter 7. 



3.7.3 Search for Gravitomagnetism 

According to GR, moving or rotating matter should produce a contribution 
to the gravitational field that is the analogue of the magnetic field of a mov- 
ing charge or a magnetic dipole. Although gravitomagnetism plays a role 
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in a variety of measured relativistic effects, it has not been seen to date, 
isolated from other post-Newtonian effects (for a discussion of the evidence 
for gravitomagnetism in solar system measurements and the binary pulsar, 
see ||99|, |100|| ). The Relativity Gyroscope Experiment (Gravity Probe B or 
GP-B) at Stanford University, in collaboration with NASA and Lockheed- 
Martin Corporation, is in the advanced stage of developing a space mission to 
detect this phenomenon directly ||71|. A set of four superconducting-niobium- 



coated, spherical quartz gyroscopes will be flown in a low polar Earth orbit, 
and the precession of the gyroscopes relative to the distant stars will be 
measured. In the PPN formalism, the predicted effect of gravitomagnetism 
is a precession (also known as the Lense-Thirring effect, or the dragging of 
inertial frames), given by 

dS/dr = f^LT X S , f^LT = -^(1 + 7 + ^ai)[J - 3n(n ■ 3)]/r^ , (44) 

where J is the angular momentum of the Earth, n is a unit radial vector, 
and r is the distance from the center of the Earth (TEGP 9.1). For a polar 
orbit at about 650 km altitude, this leads to a secular angular precession at a 
rate |(1 + 7+ \ai)42 x 10~^ arcsec/yr. The accuracy goal of the experiment 
is about 0.5 milliarcseconds per year. The science instrument package and 
the spacecraft are in the final phases of assembly, with launch scheduled for 
2002. 

Another proposal to look for an effect of gravitomagnetism is to measure 
the relative precession of the line of nodes of a pair of laser-ranged geo- 
dynamics satellites (LAGEOS), ideally with supplementary inclination an- 
gles; the inclinations must be supplementary in order to cancel the dominant 
nodal precession caused by the Earth's Newtonian gravitational multipole 
moments. Unfortunately, the two existing LAGEOS satellites are not in ap- 
propriately inclined orbits, and no plans exist at present to launch a third 
satellite in a supplementary orbit. Nevertheless, by combing nodal preces- 
sion data from LAGEOS I and II with perigee advance data from the slightly 
eccentric orbit of LAGEOS II, Ciufolini et al. reported a partial cancellation 



of multipole effects, and a resulting 20 percent confirmation of GR p3 |. 



3.7.4 Improved PPN Parameter Values 

A number of advanced space missions have been proposed in which spacecraft 
orbiters or landers and improved tracking capabilities could lead to significant 
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improvements in values of the PPN parameters, of J2 of the Sun, and of 
G/G. Doppler tracking of the Cassini spacecraft (launched to orbit and study 
Saturn in 1997) during its 2003 superior conjunction could measure 7 to a 
few parts in 10^, by measuring the time variation of the Shapiro delay 



A Mercury orbiter, in a two-year experiment, with 3 cm range capability, 
could yield improvements in the perihelion shift to a part in 10^, in 7 to 
4 X 10~^, in G/G to 10~^^ yr~^, and in J2 to a few parts in lO*^. Proposals 
are being developed, primarily in Europe, for advanced space missions which 
will have tests of PPN parameters as key components, including GAIA, a 
high-precision astrometric telescope (successor to Hipparcos), which could 
measure light-deflection and 7 to the 10~^ level [^. Nordtvedt ||103|| has 



argued that "grand fits" of large solar system range data sets, including 
ranging to Mercury, Mars and the Moon, could yield substantially improved 
measurements of PPN parameters. 



4 Strong Gravity and Gravitational Waves: 
A New Testing Ground 

4.1 Strong-field systems in general relativity 
4.1.1 Defining weak and strong gravity 

In the solar system, gravity is weak, in the sense that the Newtonian gravi- 
tational potential and related variables (?7(x, t) ~ ~ p/p ~ e ) are every- 
where much smaller than unity everywhere. This is the basis for the post- 
Newtonian expansion and for the "parametrized post-Newtonian" framework 
described in Sec. p.2| . "Strong-field" systems are those for which the simple 
IPN approximation of the PPN framework is no longer appropriate. This 
can occur in a number of situations: 

• The system may contain strongly relativistic objects, such as neu- 
tron stars or black holes, near and inside which e ~ 1, and the post- 
Newtonian approximation breaks down. Nevertheless, under some cir- 
cumstances, the orbital motion may be such that the interbody po- 
tential and orbital velocities still satisfy e <^ 1 so that a kind of post- 
Newtonian approximation for the orbital motion might work; however. 
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the strong-field internal gravity of the bodies could (especially in alter- 
native theories of gravity) leave imprints on the orbital motion. 



• The evolution of the system may be affected by the emission of gravita- 
tional radiation. The IPN approximation does not contain the effects 
of gravitational radiation back-reaction. In the expression for the met- 
ric given in Box 2, radiation back-reaction effects do not occur until 
0(e''/^) in goQ, O(e^) in g^i, and 0(e^/^) in gij. Consequently, in order 
to describe such systems, one must carry out a solution of the equations 
substantially beyond IPN order, sufficient to incorporate the leading 
radiation damping terms at 2.5PN order. 

• The system may be highly relativistic in its orbital motion, so that 
[f n-j y"^ I even for the interbody field and orbital velocity. Sys- 
tems like this include the late stage of the inspiral of binary systems of 
neutron stars or black holes, driven by gravitational radiation damp- 
ing, prior to a merger and collapse to a final stationary state. Binary 
inspiral is one of the leading candidate sources for detection by a world- 
wide network of laser interferometric gravitational-wave observatories 
nearing completion. A proper description of such systems requires not 
only equations for the motion of the binary carried to extraordinar- 
ily high PN orders (at least 3.5PN), but also requires equations for 
the far-zone gravitational waveform measured at the detector, that are 
equally accurate to high PN orders beyond the leading "quadrupole" 
approximation. 



Of course, some systems cannot be properly described by any post-Newt- 
onian approximation because their behavior is fundamentally controlled by 
strong gravity. These include the imploding cores of supernovae, the final 
merger of two compact objects, the quasinormal-mode vibrations of neutron 
stars and black holes, the structure of rapidly rotating neutron stars, and so 
on. Phenomena such as these must be analysed using different techniques. 
Chief among these is the full solution of Einstein's equations via numerical 
methods. This field of "numerical relativity" is a rapidly growing and matur- 
ing branch of gravitational physics, whose description is beyond the scope of 
this article. Another is black hole perturbation theory (see ||9^ for a review). 
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4.1.2 Compact bodies and the Strong Equivalence Principle 

When dealing with the motion and gravitational-wave generation by orbit- 
ing bodies, one finds a remarkable simplification within general relativity. As 
long as the bodies are sufficiently well-separated that one can ignore tidal 
interactions and other effects that depend upon the finite extent of the bodies 
(such as their quadrupole and higher multipole moments), then all aspects 
of their orbital behavior and gravitational wave generation can be character- 
ized by just two parameters: mass and angular momentum. Whether their 
internal structure is highly relativistic, as in black holes or neutron stars, or 
non-relativistic as in the Earth and Sun, only the mass and angular momen- 
tum are needed. Furthermore, both quantities are measurable in principle 
by examining the external gravitational field of the bodies, and make no 
reference whatsoever to their interiors. 



Damour ||3^ calls this the "effacement" of the bodies' internal structure. 
It is a consequence of the Strong Equivalence Principle (SEP), described in 
Section |3.1.2| . 



General relativity satisfies SEP because it contains one and only one 
gravitational field, the spacetime metric g^y. Consider the motion of a body 
in a binary system, whose size is small compared to the binary separation. 
Surround the body by a region that is large compared to the size of the body, 
yet small compared to the separation. Because of the general covariance of 
the theory, one can choose a freely-falling coordinate system which comoves 
with the body, whose spacetime metric takes the Minkowski form at its outer 
boundary (ignoring tidal effects generated by the companion) . There is thus 
no evidence of the presence of the companion body, and the structure of 
the chosen body can be obtained using the field equations of GR in this 
coordinate system. Far from the chosen body, the metric is characterized 
by the mass and angular momentum (assuming that one ignores quadrupole 
and higher multipole moments of the body) as measured far from the body 
using orbiting test particles and gyroscopes. These asymptotically measured 
quantities are oblivious to the body's internal structure. A black hole of mass 
m and a planet of mass m would produce identical spacetimes in this outer 
region. 

The geometry of this region surrounding the one body must be matched to 
the geometry provided by the companion body. Einstein's equations provide 
consistency conditions for this matching that yield constraints on the motion 
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of the bodies. These are the equations of motion. As a result the motion of 
two planets of mass and angular momentum mi, 1712, Ji and J2 is identical 
to that of two black holes of the same mass and angular momentum (again, 
ignoring tidal effects). 

This effacement does not occur in an alternative gravitional theory like 
scalar-tensor gravity. There, in addition to the spacetimc metric, a scalar 
field (f) is generated by the masses of the bodies, and controls the local value 
of the gravitational coupling constant {i.e. G is a function of (f). Now, in 
the local frame surrounding one of the bodies in our binary system, while 
the metric can still be made Minkowskian far away, the scalar field will 
take on a value 0o determined by the companion body. This can affect the 
value of G inside the chosen body, alter its internal structure (specifically its 
gravitational binding energy) and hence alter its mass. Effectively, each mass 
becomes several functions 771^(0) of the value of the scalar field at its location, 
and several distinct masses come into play, incrtial mass, gravitational mass, 
"radiation" mass, etc. The precise nature of the functions will depend on 
the body, specifically on its gravitational binding energy, and as a result, 
the motion and gravitational radiation may depend on the internal structure 
of each body. For compact bodies such as neutron stars, and black holes 
these internal structure effects could be large; for example, the gravitational 
binding energy of a neutron star can be 40 percent of its total mass. At IPN 
order, the leading manifestation of this effect is the Nordtvedt effect. 

This is how the study of orbiting systems containing compact objects 
provides strong-field tests of general relativity. Even though the strong-field 
nature of the bodies is effaced in GR, it is not in other theories, thus any 
result in agreement with the predictions of GR constitutes a kind of "null" 
test of strong-field gravity. 

4.2 Motion and gravitational radiation in general rel- 
ativity 

The motion of bodies and the generation of gravitational radiation are long- 
standing problems that date back to the first years following the publication 
of GR, when Einstein calculated the gravitational radiation emitted by a 
laboratory-scale object using the linearized version of GR. and de Sitter cal- 
culated N-body equations of motion for bodies in the IPN approximation 
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to GR. It has at times been controversial, with disputes over such issues 
as whether Einstein's equations alone imply equations of motion for bodies 
(Einstein, Infeld and Hoffman demonstrated explicitly that they do, using 
a matching procedure similar to the one described above), whether gravita- 
tional waves are real or are artifacts of general covariance (Einstein waffled; 
Bondi and colleagues proved their reality rigorously in the 1950s), and even 
over algebraic errors (Einstein erred by a factor of 2 in his first radiation 
calculation; Eddington found the mistake). Shortly after the discovery of the 
binary pulsar PSR 1913+16 in 1974, questions were raised about the founda- 
tions of the "quadrupole formula" for gravitational radiation damping (and 
in some quarters, even about its quantitative validity). These questions were 
answered in part by theoretical work designed to shore up the foundations of 
the quadrupole approximation, and in part (perhaps mostly) by the agree- 
ment between the predictions of the quadrupole formula and the observed 
rate of damping of the pulsar's orbit (see Section p.l| ). Damour |^ gives a 
thorough review of this subject. 

The problem of motion and radiation has received renewed interest since 
1990, with proposals for construction of large-scale laser interferometric grav- 
itational-wave observatories, such as the LIGO project in the US, VIRGO 
and GEO600 in Europe, and TAMA300 in Japan, and the realization that a 
leading candidate source of detectable waves would be the inspiral, driven by 
gravitational radiation damping, of a binary system of compact objects (neu- 
tron stars or black holes) [jl|, |127|] . The analysis of signals from such systems 
will require theoretical predictions from GR that are extremely accurate, well 
beyond the leading-order prediction of Newtonian or even post-Newtonian 
gravity for the orbits, and well beyond the leading-order formulae for gravi- 
tational waves. 

This presented a major theoretical challenge: to calculate the motion and 
radiation of systems of compact objects to very high PN order, a formidable 
algebraic task, while addressing a number of issues of principle that have 
historically plagued this subject, sufficiently well to ensure that the results 
were physically meaningful. This challenge is in the process of being met, 
so that we may soon see a remarkable convergence between observational 
data and accurate predictions of gravitational theory that could provide new, 
strong-field tests of GR. 

Here we give a brief overview of the problem of motion and gravitational 
radiation. 
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4.3 Einstein's equations in "relaxed" form 

The Einstein equations G^,^ = SvrT^j, are elegant and deceptively simple, 
showing geometry (in the form of the Einstein tensor G^;/, which is a func- 
tion of spacetime curvature) being generated by matter (in the form of the 
material stress-energy tensor T^^,. However, this is not the most useful form 
for actual calculations. For post-Newtonian calculations, a far more useful 
form is the so-called "relaxed" Einstein equations: 

□ = _i67rr"^ , (45) 

where □ = —d'^/dt^ -|- is the fiat-spacetime wave operator, /i"^ is a "grav- 
itational tensor potential" related to the deviation of the spacetime metric 
from its Minkowski form by the formula /i"^ = 77"^ — {—gY^'^g"^ ■, g is the 
determinant of gap, and a particular coordinate system has been specified 
by the deDonder or harmonic gauge condition dh"'^ /dx^ = (summation on 
repeated indices is assumed). This form of Einstein's equations bears a strik- 
ing similarity to Maxwell's equations for the vector potential A" in Lorentz 
gauge: DA" = —An J", dA'^/dx°' = 0. There is a key difference, however: 
the source on the right hand side of Eq. ( ^5D is given by the "effective" 
stress-energy pseudotensor 

r"/5 = (-^)T"^ + (167r)-^A°^ , (46) 

where A"^ is the non-linear "field" contribution given by terms quadratic 
(and higher) in /i"^ and its derivatives (see |9|], Eqs. (20.20) - (20.21) for 
formulae). In general relativity, the gravitational field itself generates gravity, 
a reflection of the nonlinearity of Einstein's equations, and in contrast to the 
linearity of Maxwell's equations. 

Equation (^5]) is exact, and depends only on the assumption that space- 
time can be covered by harmonic coordinates. It is called "relaxed" because 
it can be solved formally as a functional of source variables without specifying 
the motion of the source, in the form 



/i"^(t,x) = 4 / ^ ' ^ ''^ dV, (47) 



r"^(t - 
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where the integration is over the past fiat-spacetime null cone C of the field 
point (t, x). The motion of the source is then determined either by the 
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equation dr'^^/dx^ = (which follows from the harmonic gauge condition), 
or from the usual covariant equation of motion T°'^.f3 = 0, where the subscript 
; f3 denotes a covariant divergence. This formal solution can then be iterated 
in a slow motion {v < 1) weak-field < 1) approximation. One begins 

by substituting h^^ = into the source t"^ in Eq. (0), and solving for the 
first iterate hi^, and then repeating the procedure sufficiently many times to 
achieve a solution of the desired accuracy. For example, to obtain the IPN 
equations of motion, two iterations are needed {i.e. h'^'^ must be calculated); 
likewise, to obtain the leading gravitational waveform for a binary system, 
two iterations are needed. 

At the same time, just as in electromagnetism, the formal integral (^) 
must be handled differently, depending on whether the field point is in the far 
zone or the near zone. For field points in the far zone or radiation zone, |x| > 
A > |x'| (A is the gravitational wavelength/27r), the field can be expanded in 
inverse powers of i? = |x| in a multipole expansion, evaluated at the "retarded 
time" t — R. The leading term in is the gravitational waveform. For 
field points in the near zone or induction zone, |x| ~ |x'| < A, the field is 
expanded in powers of |x — x'| about the local time t, yielding instantaneous 
potentials that go into the equations of motion. 

However, because the source r"^ contains /i"^ itself, it is not confined 
to a compact region, but extends over all spacetime. As a result, there is a 
danger that the integrals involved in the various expansions will diverge or 
be ill-defined. This consequence of the non-linearity of Einstein's equations 
has bedeviled the subject of gravitational radiation for decades. Numerous 
approaches have been developed to try to handle this difficulty. The "post- 



Minkowskian" method of Blanchet, Damour and Iyer |]T8|, |g, ||, |T]' Hi 
solves Einstein's equations by two different techniques, one in the near zone 
and one in the far zone, and uses the method of singular asymptotic matching 
to join the solutions in an overlap region. The method provides a natural 
"regularization" technique to control potentially divergent integrals. The 
"Direct Integration of the Relaxed Einstein Equations" (DIRE) approach of 
Will, Wiseman and Pati ||151| , |1U5|| , retains Eq. ( ^7\ ) as the global solution. 



but splits the integration into one over the near zone and another over the 
far zone, and uses different integration variables to carry out the explicit 
integrals over the two zones. In the DIRE method, all integrals are finite and 
convergent. 

These methods assume from the outset that gravity is sufficiently weak 
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that < 1 and harmonic coordinates exists everywhere, including inside 

the bodies. Thus, in order to apply the results to cases where the bodies may 
be neutron stars or black holes, one relies upon the Strong Equivalence Prin- 
ciple to argue that, if tidal forces are ignored, and equations are expressed in 
terms of masses and spins, one can simply extrapolate the results unchanged 
to the situation where the bodies are ultrarelativistic. While no general proof 
of this exists, it has been shown to be valid in specific circumstances, such 
as at 2PN order in the equations of motion, and for black holes moving in a 
Newtonian background field . 



Methods such as these have resolved most of the issues that led to criti- 
cism of the foundations of gravitational radiation theory during the 1970s. 



4.4 Equations of motion and gravitational waveform 



Among the results of these approaches are formulae for the equations of 
motion and gravitational waveform of binary systems of compact objects, 
carried out to high orders in an FN expansion. Here we shall only state the 
key formulae that will be needed for this article. For example, the relative 
two-body equation of motion has the form 



m 



{- 



n 



A 



IPN 



A 



2.5PN 



+ AspM + A 



3.5PAf 



+ ...} 



where m 
n = (xi 



(48) 

mi -|- 777,2 is the total mass, r = |xi — X2I, v = Vi — V2, and 
X2)/r. The notation AnPN indicates that the term is O^e"^) 
relative to the Newtonian term — fi. Explicit formulae for terms through 
various orders have been calculated by various authors: non-radiative terms 
through 2FN order [E^, RB, pTl E6i, radiation reaction terms at 2.5FN and 



3.5FN order 78 
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and non-radiative 3FN terms ^ pq] . 
Here we quote only the first FN corrections and the leading radiation-reaction 
terms at 2.5FN order: 



AiPAT 
A2.5PAf 



|(4 + 2ri) (1 + 3r])v'^ + -7/7-H n + (4 - 27/)rv , (49) 

(50) 
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where rj = 17111712/ {mi + m2Y. These terms are sufficient to analyse the orbit 



and evolution of the binary pulsar (Sec. p.l| ). For example, the IFN terms 



59 



are responsible for the periastron advance of an eccentric orbit, given by 
uj = 67r/fem/a(l — e^), where a and e are the semi-major axis and eccentricity, 
respectively, of the orbit, and fb is the orbital frequency, given to the needed 
order by Kepler's third law 27rfb = {m/a?)^/'^. 

Another product is a formula for the gravitational field far from the sys- 
tem, written schematically in the form 

K'^ = IQ*"* + Qo,5PN + QlPN + Ql.bPN + QipN + Q2.5PN + • • •} ; (51) 

where R is the distance from the source, and the variables are to be evaluated 
at retarded time t — R. The leading term is the so-called quadrupole formula 

/i^^(t,x) = |^(t-i?), (52) 

where is the quadrupole moment of the source, and overdots denote time 
derivatives. For a binary system this leads to 

Q'^ = 27]{vV - mh'h^/r) . (53) 

For binary systems, explicit formulae for all the terms through 2.5PN order 
have been derived by various authors ||134| , |154| , |106| , p2| , |23| , |151| , p!5| , |l^]. 



Given the gravitational waveform, one can compute the rate at which energy 
is carried off by the radiation (schematically / hhdQ, the gravitational analog 
of the Poynting flux). The lowest-order quadrupole formula leads to the 
gravitational-wave energy flux 



4 

,m 



E = -^' — {12v'-llr'). (54) 

Formulae for fluxes of angular and linear momentum can also be derived. 
The 2.5PN radiation-reaction terms in the equation of motion (^) result in 
a damping of the orbital energy that precisely balances the energy flux (^) 
determined from the waveform. Averaged over one orbit, this results in a 
rate of increase of the binary's orbital frequency given by 

A = ^-^m2nMh)'/'F{e) , (55) 
5 

where A4 is the so-called "chirp" mass, given hy A4 = rf/^m, and F{e) = (1 + 
73eV24 + 37eV96) /(I - e^/^. Notice that by making precise measurements 
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of the phase = 2-k f(t')dt' of either the orbit or the gravitational 

waves (for which / = 2/^ for the dominant component) as a function of the 
frequency, one in effect measures the "chirp" mass of the system. 

These formahsms have also been generalized to include the leading effects 
of spin-orbit and spin-spin coupling between the bodies ||82|, |83| . 

Another approach to gravitational radiation is applicable to the special 
limit in which one mass is much smaller than the other. This is the method 
of black-hole perturbation theory. One begins with an exact background 
spacetime of a black hole, either the non-rotating Schwarzschild or the ro- 
tating Kerr solution, and perturbs it according to gf^u = gj^J + h^^. The 
particle moves on a geodesic of the background spacetime, and a suitably 
defined source stress-energy tensor for the particle acts as a source for the 
gravitational perturbation and wave field hfj_^. This method provides nu- 
merical results that are exact in f , as well as analytical results expressed as 
series in powers of v, both for non-rotating and for rotating black holes. For 
non- rotating holes, the analytical expansions have been carried to 5.5 PN 
order, or e^'^ beyond the quadrupole approximation. All results of black hole 
perturbation agree precisely with the mi — > limit of the PN results, up to 
the highest PN order where they can be compared (for a detailed review see 



4.5 Gravitational- wave detection 

A gravitational-wave detector can be modelled as a body of mass M at a 
distance L from a fiducial laboratory point, connected to the point by a 
spring of resonant frequency uq and quality factor Q. From the equation of 
geodesic deviation, the infinitesimal displacement ^ of the mass along the line 
of separation from its equilibrium position satisfies the equation of motion 



e + + ool^ = I (f+(^, 0, ^)h+{t) + {9, 0, ^)h, (t)) , (56) 



where F+{9,(l),ip) and F^{9,(l),'il}) are "beam-pattern" factors, that depend 
on the direction of the source {9, 0), and on a polarization angle ip, and /i+(t) 
and hx (t) are gravitational waveforms corresponding to the two polarizations 
of the gravitational wave (for a review, see |p.26|| ). In a source coordinate 



system in which the x — y plane is the plane of the sky and the 2;-direction 
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points toward the detector, these two modes are given by 

h4t) = ^ih^%t)~hyy{t)), h^{t) = h^y{t), (57) 

where h^^ represent transverse-traceless (TT) projections of the calculated 
waveform of Eq. (|5T|) , given by 

/i^'^ = h^\{5ik - N'N''){6]1 - N^N^) - ]^{5i] - N'N^){5kl - N^N^)\ , (58) 

where A^-' is a unit vector pointing toward the detector. The beam pattern 
factors depend on the orientation and nature of the detector. For a wave 
approaching along the laboratory z-direction, and for a mass whose location 
on the x — y plane makes an angle with the x axis, the beam pattern factors 
are given by = cos 20 and Fx = sin 20. For a resonant cylinder oriented 
along the laboratory z axis, and for source direction (^,0), they are given 
by = sin^^cos2'0, Fx = sin^6'sin2-0 (the angle ■0 measures the relative 
orientation of the laboratory and source x-axes). For a laser interferometer 
with one arm along the laboratory a;-axis, the other along the |/-axis, and 
with ^ defined as the differential displacement along the two arms, the beam 
pattern functions are -F4. = |(l + cos^6') cos 20 cos 2-0 — cos 6' sin 20 sin 2-0 and 
Fx = |(1 + cos^ &) cos 20 sin 2?/) + cos 6' sin 20 cos 2-0. 

The waveforms /i+(t) and hy^{^) depend on the nature and evolution of 
the source. For example, for a binary system in a circular orbit, with an 
inclination i relative to the plane of the sky, and the x-axis oriented along 
the major axis of the projected orbit, the quadrupole approximation of Eq. 
(H) gives 

K{t) = -^(27rA</,)2/3(l + cos2^)cos2$,(t), 

h^it) = -^(27rA</b)2/^2cosicos2$fc(t) , (59) 

where $fe(t) = 2n /* ff,(t')dt' is the orbital phase. 

5 Stellar system tests of gravitational theory 
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5.1 The binary pulsar and general relativity 

The 1974 discovery of the binary pulsar PSR 1913+16 by Joseph Taylor 
and Russell Hulse during a routine search for new pulsars provided the first 
possibility of probing new aspects of gravitational theory: the effects of strong 
relativistic internal gravitational fields on orbital dynamics, and the effects 
of gravitational radiation reaction. For reviews of the discovery and current 
status, see the published Nobel Prize lectures by Hulse and Taylor |76|, |123|] . 
For a thorough review of pulsars, including binary and millisecond pulsars. 



see 



The system consists of a pulsar of nominal period 59 ms in a close binary 
orbit with an as yet unseen companion. The orbital period is about 7.75 
hours, and the eccentricity is 0.617. From detailed analyses of the arrival 
times of pulses (which amounts to an integrated version of the Doppler-shift 
methods used in spectroscopic binary systems), extremely accurate orbital 
and physical parameters for the system have been obtained (Table Be- 
cause the orbit is so close (~ IRq) and because there is no evidence of an 
eclipse of the pulsar signal or of mass transfer from the companion, it is 
generally believed that the companion is compact: evolutionary arguments 
suggest that it is most likely a dead pulsar. Thus the orbital motion is very 
clean, free from tidal or other complicating effects. Furthermore, the data 
acquisition is "clean" in the sense that by exploiting the intrinsic stability of 
the pulsar clock combined with the ability to maintain and transfer atomic 
time accurately using such devices as the Global Positioning System, the 
observers can keep track of the pulsar phase with an accuracy of 15/is, de- 
spite extended gaps between observing sessions (including a several-year gap 
during the middle 1990s upgrade of the Arecibo radio telecsope) . The pulsar 
has shown no evidence of "glitches" in its pulse period. 

Three factors make this system an arena where relativistic celestial me- 
chanics must be used: the relatively large size of relativistic effects [forbit ~ 
(m/r)^/^ ~ 10~^], a factor of 10 larger than the corresponding values for 
solar-system orbits; the short orbital period, allowing secular effects to build 
up rapidly; and the cleanliness of the system, allowing accurate determina- 
tions of small effects. Because the orbital separation is large compared to 
the neutron stars' compact size, tidal effects can be ignored. Just as New- 
tonian gravity is used as a tool for measuring astrophysical parameters of 
ordinary binary systems, so GR is used as a tool for measuring astrophysical 
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Symbol 




Parameter 


(units) 


Value! 


(i) "Physical" Parameters 






Right Ascension 


a 


19''15'^28/00018(15) 


Declination 


5 


16°06'27."4043(3) 


Pulsar Period 


Pp (ms) 


59.029997929613(7) 


Derivative of Period 


Pp 


8.62713(8) X 10-18 


(ii) "Keplerian" Parameters 






Projected semimajor axis 


Op sin i (s) 


2.3417592(19) 


Eccentricity 


e 


0.6171308(4) 


Orbital Period 


Pb (day) 


0.322997462736(7) 


Longitude of periastron 


^0 n 


226.57528(6) 


Julian date of periastron 


To (MJD) 


46443.99588319(3) 



(iii) "Post-Keplerian" Parameters 

Mean rate of periastron advance {i1j){°jt^^) 4.226621(11) 

Redshift/time dilation 7' (ms) 4.295(2) 

Orbital period derivative A (lO^^^) -2.422(6) 



Numbers in parentheses denote errors in last digit. 

Data from Itittp:/ /puppsr8. princeton.edu/psrcat. html 

Table 6: Parameters of the Binary Pulsar PSR 1913+16 
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parameters in the binary pulsar. 

The observational parameters that are obtained from a least-squares solu- 
tion of the arrival-time data fall into three groups: (i) non-orbital parameters, 
such as the pulsar period and its rate of change (defined at a given epoch), 
and the position of the pulsar on the sky; (ii) five "Keplerian" parameters, 
most closely related to those appropriate for standard Newtonian binary 
systems, such as the eccentricity e and the orbital period Pi,; and (iii) five 
"post-Keplerian" parameters. The five post-Keplerian parameters are (cl;), 
the average rate of periastron advance; 7', the amplitude of delays in arrival 
of pulses caused by the varying effects of the gravitational redshift and time 
dilation as the pulsar moves in its elliptical orbit at varying distances from 
the companion and with varying speeds; Pb, the rate of change of orbital 
period, caused predominantly by gravitational radiation damping; and r and 
s = sini, respectively the "range" and "shape" of the Shapiro time delay of 
the pulsar signal as it propagates through the curved spacetime region near 
the companion, where i is the angle of inclination of the orbit relative to the 
plane of the sky. 

In GR, these post-Keplerian parameters can be related to the masses 
of the two bodies and to measured Keplerian parameters by the equations 
(TEGP 12.1, 14.6(a)) 

i = e{2'KU)-\2'Kmhf'\m2/m){l+m2/m), 

A = -(1927r/5)(27rM/;.)5/^F(e), 

s — sin i , 

r = m2 , (60) 

where m\ and m2 denote the pulsar and companion masses, respectively. The 
formula for {CS) ignores possible non-relativistic contributions to the perias- 
tron shift, such as tidally or rotationally induced effects caused by the com- 
panion (for discussion of these effects, sec TEGP 12.1(c)). The formula for 
P(, includes only quadrupole gravitational radiation; it ignores other sources 
of energy loss, such as tidal dissipation (TEGP 12.1(f)). Notice that, by 
virtue of Kepler's third law, (27r/b)^ = m/a^, (27rm/b)^/^ ~ m/a ~ e, thus 
the first two post-Keplerian parameters can be seen as 0(e), or IPN cor- 
rections to the underlying variable, while the third is an 0(6^/2), or 2.5PN 
correction. The current observed values for the Keplerian and post-Keplerian 
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parameters are shown in Table |^. The parameters r and s are not separately 
measurable with interesting accuracy for PSR 1913+16 because the orbit's 
47° inclination does not lead to a substantial Shapiro delay. 

Because ft, and e are separately measured parameters, the measurement 
of the three post-Keplerian parameters provide three constraints on the two 
unknown masses. The periastron shift measures the total mass of the sys- 
tem, Pb measures the chirp mass, and 7' measures a complicated function of 
the masses. GR passes the test if it provides a consistent solution to these 
constraints, within the measurement errors. 

From the intersection of the {u) and 7' constraints we obtain the values 
mi = 1.4411 ±0.0007Mq and ms = 1.3873 ±0.0007Mq. The third of Eqs. |60| 
then predicts the value Pb = -2.40243±0.00005 x 10^^^ Q^der to compare 
the predicted value for Pt, with the observed value of Table ||, it is necessary 
to take into account the small effect of a relative acceleration between the 
binary pulsar system and the solar system caused by the differential rotation 
of the galaxy. This effect was previously considered unimportant when Pb 



was known only to 10 percent accuracy. Damour and Taylor ||51| carried out a 
careful estimate of this effect using data on the location and proper motion of 
the pulsar, combined with the best information available on galactic rotation, 
and found 

ff^L ~ -(1.7±0.5) X 10"^^ (61) 
Subtracting this from the observed Pb (Table ^) gives the residual 

pCORR _ _(2.408 ± 0.010[OBS] ± 0.005[GAL]) x 10'^^ ^ (52) 
which agrees with the prediction, within the errors. In other words, 

pGR 

- 1.0023 ± 0.0041 [OBS] ± 0.0021 [GAL] . (63) 



pCORR 



The consistency among the measurements is displayed in Figure ^ in which 
the regions allowed by the three most precise constraints have a single com- 
mon overlap. 

A third way to display the agreement with general relativity is by com- 
paring the observed phase of the orbit with a theoretical template phase 
as a function of time. If varies slowly in time, then to first order in 
a Taylor expansion, the orbital phase is given by $b(t) = ^nfbot + Tcfbot"^- 
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The time of periastron passage tp is given by ^(tp) = 2ttN, where N 
is an integer, and consequently, the periastron time will not grow linearly 
with A^. Thus the cumulative difference between periastron time tp and 
N/ fbQ, the quantities actually measured in practice, should vary according 
to tp - N/fho = -AoArV2/3) ^ -(Ao/2/w)t'. Figure shows the results: 
the dots are the data points, while the curve is the predicted difference using 



the measured masses and the quadrupole formula for /^o ||136 . 

The consistency among the constraints provides a test of the assumption 
that the two bodies behave as "point" masses, without complicated tidal 
effects, obeying the general relativistic equations of motion including grav- 
itational radiation. It is also a test of strong gravity, in that the highly 
relativistic internal structure of the neutron stars does not influence their 
orbital motion, as predicted by the Strong Equivalence Principle of GR. 



Recent observations |137[| indicate variations in the pulse profile, which 



suggests that the pulsar is undergoing precession as it moves through the 
curved spacetime generated by its companion, an effect known as geodetic 
precession. The amount is consistent with GR, assuming that the pulsar's 
spin is suitably misaligned with the orbital angular momentum. Unfortu- 
nately, the evidence suggests that the pulsar beam may precess out of our 
line of sight by 2020. 

5.2 A population of binary pulsars? 

Since 1990, several new massive binary pulsars similar to PSR 1913+16 were 
discovered, leading to the possibility of new or improved tests of GR. 

PSR 1534+12. This is a binary pulsar system in our galaxy. Its pulses 
are significantly stronger and narrower than those of PSR 1913+16, so timing 
measurements are more precise, reaching 3/is accuracy. Its parameters are 
listed in Table |^ ||118| , |119|] . The orbital plane appears to be almost edge- 



on relative to the line of sight {i ~ 80°); as a result the Shapiro delay is 
substantial, and separate values of the parameters r and s have been obtained 
with interesting accuracy. Assuming general relativity, one infers that the 
two masses are mi = 1.335 ± O.OO2M0 and ma = 1.344 ± O.OO2M0. The rate 
of orbit decay Pf, agrees with GR to about 15 percent, the precision limited 
by the poorly known distance to the pulsar, which introduces a significant 
uncertainty into the subtraction of galactic acceleration. Independently of 
Pf), measurement of the four other post-Keplerian parameters gives two tests 
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1.42 1.43 1.44 1.45 1.46 
MASS OF PULSAR (solar masses)- 



Figure 6: Constraints on masses of pulsar and companion from data on 
PSR 1913+16, assuming GR to be valid. Width of each strip in the plane 
reflects observational accuracy, shown as a percentage. Inset shows the three 
constraints on the full mass plane; intersection region (a) has been magnified 
400 times for the full figure. 
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Figure 7: Plot of cumulative shift of the periastron time from 1975 - 2000. 
Points are data, curve is the GR prediction. Gap during the middle 1990s 
was caused by closure of Arecibo for upgrading. [J. H. Taylor and J. M. 
Weisberg, 2000, private communication]. 
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Parameter B1534+12 B2127+11C B1855+09 B0655+64 

(i) "Keplerian" Parameters 

apsini (s) 3.729464(3) 2.520(3) 9.2307802(4) 4.125612(5) 

e 0.2736775(5) 0.68141(2) 0.00002168(5) 0.0000075(11) 

Pb(day) 0.42073729933(3) 0.335282052(6) 12.3271711905(6) 1.028669703(1) 

(ii) "Post-Keplerian" Parameters^ 

((j)(°yr-i) 1.755794(19) 4.457(12) 
y (ms) 2.071(6) 4.9(1.1) 

A(10"^2) -0.131(9) < 0.5 

r(/xs) 6.3(1.3) 1.27(10) 

s = smi 0.983(8) 0.9992(5) 



1 From P§^mjm 



and [http:/ /puppsr8.princeton.edu/psrcat.htmi| 



Table 7: Parameters of Other Binary Pulsars 



of strong- field gravity in the non-radiative regime [124 



PSR 2127+1 IC. This system appears to be a clone of the Hulse- Taylor 
binary pulsar, with very similar values for orbital period and eccentricity (see 
Table |^. The inferred total mass of the system is 2.706 iO.OllM©. Because 
the system is in the globular cluster M15 (NGC 7078), it suffers Doppler shifts 
resulting from local accelerations, either by the mean cluster gravitational 
field or by nearby stars, that are more difficult to estimate than was the case 
with the galactic system PSR 1913+16. This may make a separate, precision 
measurement of the relativistic contribution to Pb impossible. 

PSR 1855+09. This binary pulsar system is not particularly relativis- 
tic, with a long period (12 days) and highly circular orbit. However, because 
we observe the orbit nearly edge on, the Shapiro delay is large and measur- 
able, as reflected in the post-Keplerian parameters r and s. 

PSR 0655+64- This system consists of a pulsar and a white dwarf 
companion in a nearly circular orbit. Only an upper limit on Ph has been 
placed. 
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5.3 Binary pulsars and alternative theories 

Soon after the discovery of the binary pulsar it was widely hailed as a new 
testing ground for relativistic gravitational effects. As we have seen in the 
case of GR, in most respects, the system has lived up to, indeed exceeded, 
the early expectations. 

In another respect, however, the system has only partially lived up to its 
promise, namely as a direct testing ground for alternative theories of grav- 
ity. The origin of this promise was the discovery that alternative theories 
of gravity generically predict the emission of dipole gravitational radiation 
from binary star systems. In general relativity, there is no dipole radiation 
because the "dipole moment" (center of mass) of isolated systems is uniform 
in time (conservation of momentum), and because the "inertial mass" that 
determines the dipole moment is the same as the mass that generates grav- 
itational waves (SEP). In other theories, while the inertial dipole moment 
may remain uniform, the "gravity-wave" dipole moment need not, because 
the mass that generates gravitational waves depends differently on the in- 
ternal gravitational binding energy of each body than does the inertial mass 
(violation of SEP). Schematically, in a coordinate system in which the center 
of inertial mass is at the origin, so that m/^iXi -|- mi ^2^2 = 0, the dipole part 
of the retarded gravitational field would be given by 



where v = Vi — V2 and rj and m are defined using inertial masses. In theories 
that violate SEP, the difference between gravitational-wave mass and inertial 
mass is a function of the internal gravitational binding energy of the bod- 
ies. This additional form of gravitational radiation damping could, at least 
in principle, be significantly stronger than the usual quadrupole damping, 
because it depends on fewer powers of the orbital velocity v, and it depends 
on the gravitational binding energy per unit mass of the bodies, which, for 
neutron stars, could be as large as 40 percent (see TEGP 10 for further de- 
tails). As one fulfillment of this promise. Will and Eardley worked out in 
detail the effects of dipole gravitational radiation in the bimetric theory of 
Rosen, and, when the first observation of the decrease of the orbital period 
was announced in 1979, the Rosen theory suffered a terminal blow. A wide 
class of alternative theories also fail the binary pulsar test because of dipole 




(64) 
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gravitational radiation (TEGP 12.3). 

On the other hand, the early observations of PSR 1913+16 already indi- 
cated that, in GR, the masses of the two bodies were nearly equal, so that, in 
theories of gravity that are in some sense "close" to GR, dipole gravitational 
radiation would not be a strong effect, because of the apparent symmetry of 
the system. The Rosen theory, and others like it, are not "close" to general 
relativity, except in their predictions for the weak-field, slow-motion regime 
of the solar system. When relativistic neutron stars are present, theories 
like these can predict strong effects on the motion of the bodies resulting 
from their internal highly relativistic gravitational structure (violations of 
SEP). As a consequence, the masses inferred from observations of the peri- 
astron shift and 7' may be significantly different from those inferred using 
general relativity, and may be different from each other, leading to strong 
dipole gravitational radiation damping. By contrast, the Brans-Dicke theory 
is "close" to GR, roughly speaking within I/cubd of the predictions of the 
latter, for large values of the coupling constant cjbd (here we use the sub- 
script BD to distinguish the coupling constant from the periastron advance 
u). Thus, despite the presence of dipole gravitational radiation, the binary 
pulsar provides at present only a weak test of Brans-Dicke theory, not yet 
competitive with solar-system tests. 

5.4 Binary pulsars and scalar-tensor gravity 

Making the usual assumption that both members of the system are neutron 
stars, and using the methods summarized in TEGP Chapters 10-12, one can 
obtain formulas for the periastron shift, the gravitational redshift /second- 
order Doppler shift parameter, and the rate of change of orbital period, anal- 
ogous to Eqs. (^). These formulas depend on the masses of the two neutron 
stars, on their self-gravitational binding energy, represented by "sensitivities" 
s and and on the Brans-Dicke coupling constant ubd- First, there is a 
modification of Kepler's third law, given by 

2nfb = {gm/a'y/\ (65) 

Then, the predictions for (uj), 7' and Pb are 

{dj) = QTxh{2T:mhf'\l-e^r^Vg-^''' , (66) 
7' = e{2'Kh)-\2Timhfl%m2/m)g-^'^{al + gm2/m + K[r,l), (67) 
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A = -(1927r/5)(27rA^/,)5/3F'(e)-47r(27r/x/,)e52G(e), (68) 



where A4 = x^^^G ^l^rfl^m^ and, to first order in ^ = (2 + cubd) we have 



F\e) = 






G{e) = 


(l-eV/^(l + ie^), 




S = 


Sl - S2 , 




Q = 


1 - ^{Sl + S2- 2S1S2) , 




V = 




-2S1S2)], 


* 

"2 = 


1-^S2, V*2 = (1 - 


2s2)e, 


X = 






r = 


1 — 2(miS2 + m2Si)/m 


, r' = 1 - si - S2 



The quantities Sa and are defined by 

/9(lnm„)\ ^ /9(ln4 



(69) 
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and measure the "sensitivity" of the mass and moment of inertia of 
each body to changes in the scalar field (refiected in changes in G) for a 
fixed baryon number (see TEGP 11, 12 and 14.6(c) for further details). 
The quantity Sa is related to the gravitational binding energy. Notice how 
the violation of SEP in Brans-Dicke theory introduces complex structure- 
dependent effects in everything from the Newtonian limit (modification of the 
effective coupling constant in Kepler's third law) to gravitational-radiation. 
In the limit ,^ — > 0, we recover GR, and all structure dependence disappears. 
The first term in [Eq. (|5BD] is the effect of quadrupole and monopole 
gravitational radiation, while the second term is the effect of dipole radiation. 

In order to estimate the sensitivities Sa and k*, one must adopt an equa- 
tion of state for the neutron stars. It is sufficient to restrict attention to 
relatively stiff neutron star equations of state in order to guarantee neutron 
stars of sufficient mass, approximately 1.4M0. The lower limit on uj-qd re- 
quired to give consistency among the constraints on (a;), 7 and Pb as in Figure 



is several hundred [p,52||. The combination of {Cj) and 7 give a constraint on 
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the masses that is relatively weakly dependent on ^, thus the constraint on ^ 
is dominated by Pb and is directly proportional to the measurement error in 
Pb] in order to achieve a constraint comparable to the solar system value of 
3 X 10~^, the error in P^^^ would have to be reduced by more than a factor 
of ten. 

Alternatively, a binary pulsar system with dissimilar objects, such as a 
white dwarf or black hole companion, would provide potentially more promis- 
ing tests of dipole radiation. Unfortunately, none has been discovered to date; 
the dissimilar system B0655+64, with a white dwarf companion is in a highly 
circular orbit, making measurement of the periastron shift meaningless, and 
is not as relativistic as 1913+16. From the upper limit on Pb (Table one 
can infer at best the weak bound ljbd > 100 



Damour and Esposito-Farese [|T| have generalized these results to a broad 
class of scalar-tensor theories. These theories are characterized by a single 
function a{ip) of the scalar field (p, which mediates the coupling strength of 
the scalar field. For application to the solar system or to binary systems, one 
expands this function about a cosmo logical background field value ipo: 

= ao{f - ipo) + ^Poif - ipoY + ... . (71) 

A purely linear coupling function produces Brans-Dicke theory, with Oq = 
1/{2ubd + 3). The function a{ip) acts as a potential function for the scalar 
field (f, and, if /?o > 0, during cosmological evolution, the scalar field nat- 
urally evolves toward the minimum of the potential, i.e. toward ~ 0, 



ujbd oo, or toward a theory close to, though not precisely GR ^ 



Bounds on the parameters ao and /3o from solar-system, binary-pulsar and 



gravitational- wave observations (see Sec. |6.3|) are shown in Figure |g [^]. 
Negative values of Pq correspond to an unstable scalar potential; in this case, 
objects such as neutron stars can experience a "spontaneous scalarization" , 
whereby the interior values of if can take on values very different from the 
exterior values, through non-linear interactions between strong gravity and 
the scalar field, dramatically affecting the stars' internal structure and the 
consequent violations of SEP. On the other hand, /3o < is of little practical 
interest, because, with an unstable ip potential, cosmological evolution would 
presumably drive the system away from the peak where ao ~ 0, toward pa- 
rameter values that could easily be excluded by solar system experiments. 
On the ao — Po plane shown in Figure p|, the ao axis corresponds to pure 
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Figure 8: Region of the scalar-tensor theory ao — (3o plane allowed by solar- 
system, binary-pulsar, and future gravitational-wave observations. A poly- 
tropic equation of state for the neutron stars was assumed. The shaded region 
is that allowed by all tests. For positive values of Po, solar-system bounds 
(labelled IPN) still are the best. (From [Damour and Esposito-Farese 1998], 
(c)1998 by the American Physical Society, reproduced by permission) 

Brans-Dicke theory, while the origin corresponds to pure GR. As discussed 
above, solar system bounds (labelled "IPN" in Figure]^) still beat the binary 
pulsars. The bounds labelled "LIGO- VIRGO" are discussed in Sec. 

6 Gravitational- wave tests of gravitational the- 
ory 

6.1 Gravitational- wave observatories 

Some time in the next decade, a new opportunity for testing relativistic grav- 
ity will be realized, with the commissioning and operation of kilometer-scale, 
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laser interferometric gravitational-wave observatories in the U.S. (LIGO pro- 
ject), Europe (VIRGO and GEO600 projects) and Japan (TAMA300 project). 
Gravitational-wave searches at these observatories are scheduled to com- 
mence around 2002. The LIGO broad-band antennae will have the capability 
of detecting and measuring the gravitational waveforms from astronomical 
sources in a frequency band between about 10 Hz (the seismic noise cutoff) 
and 500 Hz (the photon counting noise cutoff), with a maximum sensitivity 
to strain at around 100 Hz of /i ~ A/// ~ 10~^^ (rms). The most promising 
source for detection and study of the gravitational-wave signal is the "inspi- 
ralling compact binary" - a binary system of neutron stars or black holes 
(or one of each) in the final minutes of a death dance leading to a violent 
merger. Such is the fate, for example, of the Hulse- Taylor binary pulsar 
PSR 1913+16 in about 300 million years. Given the expected sensitivity of 
the "advanced LIGO" (around 2007), which could see such sources out to 
hundreds of megaparsecs, it has been estimated that from 3 to 100 annual 
inspiral events could be detectable. Other sources, such as supernova core 
collapse events, instabilities in rapidly rotating nascent neutron stars, signals 
from non-axisymmetric pulsars, and a stochastic background of waves, may 
be detectable (for reviews, see [|l], |127|| ; for updates on the status of various 
projects, see ||6^, |3T[|). 

A similar network of cryogenic resonant-mass gravitational antennae have 
been in operation for many years, albeit at lower levels of sensitivity {h ~ 
10"^^). While modest improvements in sensitivity may be expected in the 
future, these resonant detectors are not expected to be competitive with 
the large interferometers, unless new designs involving bars of spherical, or 
nearly spherical shape come to fruition. These systems are primarily sensitive 
to waves in relatively narrow bands about frequencies in the hundreds to 
thousands of Hz range [|104|, ^ 0, |lTg]. 

In addition, plans are being developed for an orbiting laser interferometer 
space antenna (LISA for short). Such a system, consisting of three spacecraft 
separated by millions of kilometers, would be sensitive primarily in the very 
low frequency band between 10~^ and 10~^ Hz, with peak strain sensitivity 
of order h ~ lO'^s 

In addition to opening a new astronomical window, the detailed obser- 
vation of gravitational waves by such observatories may provide the means 
to test general relativistic predictions for the polarization and speed of the 
waves, and for gravitational radiation damping. 



76 



6.2 Polarization of gravitational waves 

A laser-interferometric or resonant bar gravitational- wave detector measures 
the local components of a symmetric 3x3 tensor which is composed of 
the "electric" components of the Riemann curvature tensor, Roioj, via the 
equation of geodesic deviation, given, for a pair of freely falling particles by 
= —RoiojX^ , where denotes the spatial separation. In general there are 
six independent components, which can be expressed in terms of polarizations 
(modes with specific transformation properties under rotations and boosts). 
Three are transverse to the direction of propagation, with two represent- 
ing quadrupolar deformations and one representing a monopole "breathing" 
deformation. Three modes are longitudinal, with one an axially symmetric 
stretching mode in the propagation direction, and one quadrupolar mode 
in each of the two orthogonal planes containing the propagation direction. 
Figure ^ shows the displacements induced on a ring of freely falling test par- 
ticles by each of these modes. General relativity predicts only the first two 
transverse quadrupolar modes (a) and (b) independently of the source; these 
correspond to the waveforms h+ and hx discussed earlier (note the cos2(/) and 
sin 2(j) dependences of the displacements) . Scalar-tensor gravitational waves 
can in addition contain the transverse breathing mode (c). More general met- 
ric theories predict additional longitudinal modes, up to the full complement 
of six (TEGP 10.2). 

A suitable array of gravitational antennas could delineate or limit the 
number of modes present in a given wave. The strategy depends on whether 
or not the source direction is known. In general there are eight unknowns (six 
polarizations and two direction cosines), but only six measurables (-RoiOj)- If 
the direction can be established by either association of the waves with op- 
tical or other observations, or by time-of-flight measurements between sep- 
arated detectors, then six suitably oriented detectors suffice to determine 
all six components. If the direction cannot be established, then the system 
is underdetermined, and no unique solution can be found. However, if one 
assumes that only transverse waves are present, then there are only three un- 
knowns if the source direction is known, or five unknowns otherwise. Then 
the corresponding number (three or five) of detectors can determine the po- 
larization. If distinct evidence were found of any mode other than the two 
transverse quadrupolar modes of GR, the result would be disastrous for GR. 
On the other hand, the absence of a breathing mode would not necessarily 
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rule out scalar-tensor gravity, because the strength of that mode depends on 
the nature of the source. 

Some of the details of implementing such polarization observations have 
been worked out for arrays of resonant cylindrical, disk-shaped, spherical and 
truncated icosahedral detectors (TEGP 10.2, for recent reviews see |132|| ); 
initial work has been done to assess whether the ground-based or space- 
based laser-interferometers (or combinations of the two types) could perform 
interesting polarization measurements ||133| , |32| , |90| , |67| . Unfortunately for 



this purpose, the two LIGO observatories (in Washington and Louisiana 
states, respectively) have been constructed to have their respective arms 
as parallel as possible, apart from the curvature of the Earth; while this 
maximizes the joint sensitivity of the two detectors to gravitational waves, it 
minimizes their ability to detect two modes of polarization. 



6.3 Gravitational radiation back-reaction 

In the binary pulsar, a test of OR was made possible by measuring at least 
three relativistic effects that depended upon only two unknown masses. The 
evolution of the orbital phase under the damping effect of gravitational ra- 
diation played a crucial role. Another situation in which measurement of 
orbital phase can lead to tests of GR is that of the inspiralling compact bi- 
nary system. The key differences are that here gravitational radiation itself 
is the detected signal, rather than radio pulses, and the phase evolution alone 
carries all the information. In the binary pulsar, the first derivative of the 
binary frequency, fb, was measured; here the full nonlinear variation of fb as 
a function of time is measured. 

Broad-band laser interferometers are especially sensitive to the phase evo- 
lution of the gravitational waves, which carry the information about the or- 
bital phase evolution. The analysis of gravitational-wave data from such 
sources will involve some form of matched filtering of the noisy detector out- 
put against an ensemble of theoretical "template" waveforms which depend 
on the intrinsic parameters of the inspiralling binary, such as the component 
masses, spins, and so on, and on its inspiral evolution. How accurate must 
a template be in order to "match" the waveform from a given source (where 
by a match we mean maximizing the cross-correlation or the signal-to-noise 
ratio)? In the total accumulated phase of the wave detected in the sensitive 
bandwidth, the template must match the signal to a fraction of a cycle. For 
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Figure 9: Six polarization modes for gravitational waves permitted in any 
metric theory of gravity. Shown is the displacement that each mode induces 
on a ring of test particles. The wave propagates in the +z direction. There 
is no displacement out of the plane of the picture. In (a), (b) and (c), the 
wave propagates out of the plane; in (d), (e), and (f), the wave propagates in 
the plane. In general relativity, only (a) and (b) are present; in scalar-tensor 
gravity, (c) may also be present. 
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two inspiralling neutron stars, around 16,000 cycles should be detected dur- 
ing the final few minutes of inspiral; this implies a phasing accuracy of 10~^ 
or better. Since f ~ 1/10 during the late inspiral, this means that correction 
terms in the phasing at the level of or higher are needed. More formal 



analyses confirm this intuition |3J, |59|, |108 



Because it is a slow-motion system {v ~ 10~ ), the binary pulsar is sen- 
sitive only to the lowest-order effects of gravitational radiation as predicted 
by the quadrupole formula. Nevertheless, the first correction terms of order 
V and to the quadrupole formula, were calculated as early as 1976 ( [|134|] , 
see TEGP 10.3). 

But for laser-interferometric observations of gravitational waves, the bot- 
tom line is that, in order to measure the astrophysical parameters of the 
source and to test the properties of the gravitational waves, it is necessary to 
derive the gravitational waveform and the resulting radiation back-reaction 
on the orbit phasing at least to 2PN order beyond the quadrupole approxi- 
mation, and probably to 3PN order. 

For the special case of non-spinning bodies moving on quasi-circular orbits 
{i.e. circular apart from a slow inspiral), the evolution of the gravitational 
wave frequency / = 2/^ through 2PN order has the form 



/ = ^f(vr>l/)^/^[l-fS? + ^^Vw)^/^ + 47r(vrm/) 
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(72) 



where t] = mim2/m'^- The first term is the quadrupole contribution [compare 
Eq. ([55|)], the second term is the IPN contribution, the third term, with the 
coefficient An, is the "tail" contribution, and the fourth term is the 2PN con- 



tribution, first reported jointly by Blanchet et al. |2^, |151|| . Calculation 
of the higher-order contributions is nearing completion. 

Similar expressions can be derived for the loss of angular momentum and 
linear momentum. (For explicit formulas for non-circular orbits, see |[70|| ) . 
These losses react back on the orbit to circularize it and cause it to inspiral. 
The result is that the orbital phase (and consequently the gravitational-wave 
phase) evolves non-linearly with time. It is the sensitivity of the broad-band 
LIGO and VIRGO-type detectors to phase that makes the higher-order con- 
tributions to df/dt so observationally relevant. A ready-to-use set of formulae 



80 



for the 2PN gravitational waveform template, including the non-linear evo- 
lution of the gravitational-wave frequency (not including spin effects) have 
been published and incorporated into the Gravitational Radiation Anal- 



ysis and Simulation Package (GRASP), a software toolkit used in LIGO. 

If the coefficients of each of the powers of / in Eq. ([72|) can be measured, 
then one again obtains more than two constraints on the two unknowns 
mi and 1712, leading to the possibility to test GR. For example, Blanchet 
and Sathyaprakash p8|, ^ have shown that, by observing a source with a 



sufficiently strong signal, an interesting test of the Air coefficient of the "tail" 
term could be performed. 

Another possibility involves gravitational waves from a small mass orbit- 
ing and inspiralling into a (possibly supermassive) spinning black hole. A 
general non-circular, non-equatorial orbit will precess around the hole, both 
in periastron and in orbital plane, leading to a complex gravitational wave- 
form that carries information about the non-spherical, strong-field spacetime 
around the hole. According to GR, this spacetime must be the Kerr space- 
time of a rotating black hole, uniquely specified by its mass and angular 
momentum, and consequently, observation of the waves could test this fun- 



damental hypothesis of GR ||1 14|, |107 



Thirdly, the dipole gravitational radiation predicted by scalar-tensor the- 
ories will result in a modification of the gravitational-radiation back-reaction, 
and thereby of the phase evolution. Including only the leading quadrupole 
and dipole contributions, one obtains, in Brans-Dicke theory. 



(73) 



where A4 = {x^^^Q~'^^^)r]'^^^m, and b is the coefficient of the dipole term, 
given hj b = (5/48)(x~^^^''^)^'5^, where Q, <S are given by Eqs. (|69D, and 
= 1/(2 + ujbd)- The effects are strongest for systems involving a neutron 
star and a black hole. Double neutron star systems are less promising because 
the small range of masses available near 1.4 Mq results in suppression of 
dipole radiation by symmetry (the sensitivity s turns out to be a relatively 
weak function of mass near 1.4 M©, for typical equations of state). For black 
holes, s = 0.5 identically, consequently double black- hole systems turn out 
to be observationally identical in the two theories. 

But for a 1.4 Mq neutron star and a 10 Mq (3 Mq) black hole at 200 Mpc, 
the bound on c^bd could be 600 (1800) (using advanced LIGO noise curves). 
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The bound increases linearly with signal-to-noise ratio [p.46|| . If one demands 



that this test be performed annually, thus requiring observation of frequent, 
and therefore more distant, weaker sources, the bounds on cjbd will be too 
weak to compete with existing solar-system bounds (this corresponds to the 
"LIGO- VIRGO" bound on the cto axis in Figure ||, which assumes a signal-to- 
noise ratio of 10). However, if one is prepared to wait 10 years for the lucky 
observation of a nearby, strong source, the resulting bound could exceed the 
current solar-system bound. The bounds are illustrated in Figure |1^ by the 
curves marked = 1 and = 1/10. Figure |10| assumes a double- neutron- 
star inspiral rate of 10~^ per year per galaxy, and a black-hole-neutron-star 
rate (3 times that, where j3 is highly uncertain. Values of /? = 1/10 , 1, and 10 
are shown. For the general class of scalar-tensor theories, the corresponding 
bounds are plotted on the — /?o plane in Figure ||, under the restrictive 
assumption of a signal-to-noise ratio (S/N) of 10 P3[. All other factors being 



equal, the bound achievable on the cto— /?o parameters of scalar-tensor gravity 
is inversely proportional to S/N. 

6.4 Speed of gravitational waves 

According to GR, in the limit in which the wavelength of gravitational waves 
is small compared to the radius of curvature of the background spacetime, the 
waves propagate along null geodesies of the background spacetime, i.e. they 
have the same speed, c, as light (in this section, we do not set c = 1). In other 
theories, the speed could differ from c because of coupling of gravitation to 
"background" gravitational fields. For example, in the Rosen bimetric theory 
with a flat background metric 77, gravitational waves follow null geodesies of 
77, while light follows null geodesies of g (TEGP 10.1). 

Another way in which the speed of gravitational waves could differ from 
c is if gravitation were propagated by a massive field (a massive graviton), 
in which case, Vg would be given by, in a local inertial frame, 

2 2 4 

3 = l-^> (74) 



c 



where irig and E are the graviton rest mass and energy, respectively. For a 
recent review of the idea of a massive graviton along with a model theory, 
see [p. 3 11 . 
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BLACK HOLE MASS 

Figure 10: Bounds on the scalar-tensor coupling constant cubd from 
gravitational- wave observations of inspiralling black-hole neutron-star sys- 
tems. The solar system bound is around wbd = 3500. 
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The most obvious way to test this is to compare the arrival times of a 
gravitational wave and an electromagnetic wave from the same event, e.g. a 
supernova. For a source at a distance the resulting value of the difference 

1 — Vg/c is 

where At = At a — (1 + Z)Ate is the "time difference", where At a and Ate 
are the differences in arrival time and emission time, respectively, of the two 
signals, and Z is the redshift of the source. In many cases. Ate is unknown, 
so that the best one can do is employ an upper bound on Atg based on 
observation or modelling. The result will then be a bound on 1 — Vg/c. 

For a massive graviton, if the frequency of the gravitational waves is such 
that hf ^ rUgC^, where h is Planck's constant, then Vg/c ^ 1 — ^{c/XgfY, 
where = h/rrigC is the graviton Compton wavelength, and the bound on 
1 — Vg/c can be converted to a bound on Xg, given by 

,2. f D 100Hz\'/V 1 \'/' 

A„ > 3 X 10^2 km — — - . (76 

^ \^200 Mpc / ; \fAtJ ^ ^ 

The foregoing discussion assumes that the source emits both gravitational 
and electromagnetic radiation in detectable amounts, and that the relative 
time of emission can be established to sufficient accuracy, or can be shown 
to be sufficiently small. 

However, there is a situation in which a bound on the graviton mass 
can be set using gravitational radiation alone [|148|| . That is the case of 
the inspiralling compact binary. Because the frequency of the gravitational 
radiation sweeps from low frequency at the initial moment of observation to 
higher frequency at the final moment, the speed of the gravitons emitted will 
vary, from lower speeds initially to higher speeds (closer to c) at the end. This 
will cause a distortion of the observed phasing of the waves and result in a 
shorter than expected overall time At^ of passage of a given number of cycles. 
Furthermore, through the technique of matched filtering, the parameters of 
the compact binary can be measured accurately, (assuming that GR is a good 
approximation to the orbital evolution, even in the presence of a massive 
graviton), and thereby the emission time Ate can be determined accurately. 
Roughly speaking, the "phase interval" /At in Eq. (^) can be measured to 
an accuracy 1/p, where p is the signal-to- noise ratio. 
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Thus one can estimate the bounds on Xg achievable for various compact 
inspiral systems, and for various detectors. For stellar-mass inspiral (neutron 
stars or black holes) observed by the LIGO/VIRGO class of ground-based 
interferometers, D ^ 200Mpc, / ^ lOOHz, and /At ~ p-^ ^ 1/10. The 
result is > lO^^km. For supermassive binary black holes (10^ to IO^Mq) 
observed by the proposed laser-interferometer space antenna (LISA), D ^ 
3Gpc, / ^ lO-^Hz, and fAt ~ p'^ ^ 1/1000. The result is Xg > 10^^ km. 

A full noise analysis using proposed noise curves for the advanced LIGO 
and for LISA weakens these crude bounds by factors between two and 10. 
These potential bounds can be compared with the solid bound Xg > 2.8 x 
10^^ km, ||122|| derived from solar system dynamics, which limit the presence 
of a Yukawa modification of Newtonian gravity of the form 

V{r) = (GM/r) exp(-r/Ag) , (77) 

and with the model-dependent bound Ag > 6 x 10^^ km from consideration 
of galactic and cluster dynamics ||131 . 



6.5 Other strong-gravity tests 

One of the central difficulties of testing general relativity in the strong-field 
regime is the possibility of contamination by uncertain or complex physics. In 
the solar system, weak-field gravitational effects could in most cases be mea- 
sured cleanly and separately from non-gravitational effects. The remarkable 
cleanliness of the binary pulsar permitted precise measurements of gravita- 
tional phenomena in a strong-field context. 

Unfortunately, nature is rarely so kind. Still, under suitable conditions, 
qualitative and even quantitive strong-field tests of general relativity can be 
carried out. 

One example is in cosmology. From a few seconds after the big bang 
until the present, the underlying physics of the universe is well understood, 
although significant uncertainties remain (amount of dark matter, value of 
the cosmological constant, the number of light neutrino families, etc. .). Some 
alternative theories of gravity that are qualitatively different from GR fail to 
produce cosmologies that meet even the minimum requirements of agreeing 
qualitatively with big-bang nucleosynthesis (BBN) or the properties of the 
cosmic microwave background (TEGP 13.2). Others, such as Brans-Dicke 
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theory, are sufficiently close to GR (for large enough ubd) that they conform 
to all cosmological observations, given the underlying uncertainties. The gen- 
eralized scalar-tensor theories, however, could have small ujbd early times, 
while evolving through the attractor mechanism to large ujbd today. One 
way to test such theories is through big-bang nucleosynthesis, since the abun- 
dances of the light elements produced when the temperature of the universe 
was about 1 MeV are sensitive to the rate of expansion at that epoch, which 
in turn depends on the strength of interaction between geometry and the 
scalar field. Because the universe is radiation-dominated at that epoch, un- 
certainties in the amount of cold dark matter or of the cosmological constant 
are unimportant. The nuclear reaction rates are reasonably well understood 
from laboratory experiments and theory, and the number of light neutrino 
families (3) conforms to evidence from particle accelerators. Thus, within 
modest uncertainties, one can assess the quantitive difference between the 
BBN predictions of GR and scalar-tensor gravity under strong-field condi- 
tions and compare with observations. The most sophisticated recent analysis 
[0 places bounds on the parameters and [3q of the generalized framework 
of Damour and Esposito-Farese (see Sec. |5.4| and Fig. H) that are weaker than 
solar-system bounds for (3q < 0.3, but substantially stronger for [3q > 0.3. 

Another example is the exploration of the spacetime near black holes 
via accreting matter. Observations of low-luminosity binary X-ray sources 
suggest that a form of accretion known as advection-dominated accretion 
flow (ADAF) may be important. In this kind of flow, the accreting gas is 
too thin to radiate its energy efficiently, but instead transports (advects) it 
inward toward the central object. If the central object is a neutron star, 
the matter hits the surface and radiates the energy away; if it is a black 
hole, the matter and its advected energy disappear. Systems in which the 
accreting object is believed to be a black hole from estimates of its mass 
are indeed observed to be under luminous, compared to systems where the 
object is believe to be a neutron star. This has been regarded as the ffist 
astrophysical evidence for the existence of black hole event horizons (for a 
While supporting one of the critical strong-field predictions 



review, see 



of GR, the observations and models are not likely any time soon to be able 
to distinguish one gravitational theory from another (except for theories that 
do not predict black holes at all). 

Another example involving accretion purports to explore the strong-field 
region just outside massive black holes in active galactic nuclei. Here, iron 
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in the inner region of a thin accretion disk is irradiated by X-ray-emitting 
material above or below the disk, and fluoresces in the Ka line. The spectral 
shape of the line depends on relativistic Doppler and curved-spacetime effects 
as the iron orbits the black hole near the innermost stable circular orbit, and 
could be used to determine whether the hole is a non-rotating Schwarzschild 
black hole, or a rotating Kerr black hole. Because of uncertainties in the 
detailed models, the results are inconclusive to date, but the combination 
of higher-resolution observations and better modelling could lead to striking 
tests of strong-field predictions of GR. 

7 Conclusions 

We find that general relativity has held up under extensive experimental 
scrutiny. The question then arises, why bother to continue to test it? One 
reason is that gravity is a fundamental interaction of nature, and as such 
requires the most solid empirical underpinning we can provide. Another is 
that all attempts to quantize gravity and to unify it with the other forces 
suggest that the standard general relativity of Einstein is not likely to be the 
last word. Furthermore, the predictions of general relativity are fixed; the 
theory contains no adjustable constants so nothing can be changed. Thus 
every test of the theory is either a potentially deadly test or a possible probe 
for new physics. Although it is remarkable that this theory, born 80 years ago 
out of almost pure thought, has managed to survive every test, the possibility 
of finding a discrepancy will continue to drive experiments for years to come. 
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